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ABSTRACT 


The Korean Air Force has determined that repairables management is one of the 
areas in which attention could be expected to lead to substantial improvement in the 
efficient management of defense resources and in maintaining an adequate level of force 
effectiveness. This thesis reviews various inventory models for the management of re- 
pairable items and develops an inventory model for the Korean Air Force. It discusses 
the characteristics of each model, and, identifies and explains the differences in each 
model with respect to assumptions, objectives, constraints, and optimization methods. 


Also this research evaluates the proposed model using sample data sets. 
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1, INTRODUCTION 
A. BACKGROUND 


The objective of the Korean Air Force supply system is to ensure the desired level 
of operational availability for its population of repairable equipment and maintain suf- 
ficient stock levels for replacement components and repair parts to support the mainte- 
nance of existing weapon systems and equipment. Also, the Korean Air Force supply 
system 1s responsible for stocking a number of items which are failure prone and have 
been designated as repairable. Because of engineering and economic considerations, at- 
tempts are made by Korean Air Force repair activities to restore these items to service- 
able condition whenever they fail. 

The repairable inventory system differs from that of the consumable inventory sys- 
tem in two ways. First, the repairable svstem contains two distinct inventories, one of 
which has items that are in a Ready-For-Issue (RFI) state and another which has items 
that are in a Non-Ready-For-Issue (NRFI) state. The first inventory contains items 
which are usable and the second contains items that must be repaired before they can 
be used. Second, the RFI inventory is made up of a mixture of new items and items that 
have been used, failed, repaired, and are ready to be used again. This is one reason why 
repairable inventory management is a complicated process. [Ref. 1: p. 2] 

The Korean Air Force Logistics Command (KAFLC) tries to determine an appro- 
priate stock level for repairables based on a primitive mathematical model. This model 
estimates stock levels using rules of thumb based on expert opinion. However, it is very 
difficult to determine an optimal stock level for a repairable without using a more so- 
phisticated inventory management model. The result of the lack of such a model 1s that 
the Korean Air Force supply system tends to maintain relatively high stock levels for 
repairable components. Although, repairable items account for only a small percentage 
of the quantitv of items stocked throughout the Korean Air Force, they account for a 
large portion of dollars invested in inventory. The characteristics of high cost and high 
essentialitv that these items typically possess make efficient inventory control difficult 
and extremelv critical. 

Historically, most inventory models have been developed for the private sector 
where the profit motive 1s important. Such models consider the various average annual 


variable costs associated with managing inventories and strive to minimize the sum of 


these costs. The tvpical relevant costs in most inventory models are ordering cost, 
holding cost and stockout cost. Certainly, these cost parameters are very important in 
the common concept of inventory management. [Ref. 2: p. 1] 

However, in military supply systems, it is hard to estimate the costs associated with 
stockouts. In addition, the military supplv system 1s not interested in profit 
maximization or cost minimization. Instead, such an organization is usually interested 
in maximizing the abilitv of its forces to respond to any threat. An objective which 


maximizes some measure of readiness with given resources 1s therefore appropriate. 


B. PURPOSE OF RESEARCH 

The effective accomplishment of each flying mission 15 dependent upon having suf- 
ficient numbers of aircraft ready to flv and to perform at their fullest capabilitv. To 
support this goal, the Korean Air Force emphasizes an extensive system of supply and 
maintenance capabilities. The Korean Air Force also tries to maximize the availability 
of aircraft by quickly identifying malfunctioning parts, removing them, and rapidly in- 
stalling replacements which have been positioned at the base supply department. 

The purpose of this thesis is the development of an inventory model for the man- 
agement of repairable items at the wholesale level by the Korean Air Force Logistics 
Command (KAFLC). This will be accomplished by the analysis of various repairable 
inventory management models which have been developed and applied to the manage- 
ment of inventories of repairable itenis in the United States military. And also, a nu- 
merical analvsis will be accomplished based on a proposed model. Finally the inventory 


svstem parameters, and management concept will be studied. 


C PREVIEM 

Chapter II of this thesis presents an overview of the current Korean Air Force Lo- 
gistics Command (KAFLC) repairable items management process. This overview will 
explain the organization of the Korean Air Force Logistics Command, describe the 
present general methodology used for repairable items management and the current 
overall approach to the problem. 

Chapter III gives a detailed description of the specific mathematical inventory 
models in use today for the management of repairables. Some explanation of their 
mathematical approach in solving the inventory problem is given. 

In Chapter IV, the development of a model for the Korean Air Force repairables 
inventory management system is described. Model formulation and solution techniques 


are provided. 


Chapter V will show numerical results from a computer program for comparing the 
proposed model and describe major components based on the findings in Chapters III 
and IV. 

Finally, Chapter VI contains a summary and recommendations for future studv 


based on the proposed repairables inventory management model. 


П. REPAIRABLES AND THE KOREAN AIR FORCE LOGISTICS 
SYSTEM 


A. INTRODUCTION 

The Korean Air Force didn’t use any kind of repairables inventory management 
model until the late 1960’s. In the beginning of the 1970's, the first inventory manage- 
ment concept was introduced by the U.S Air Force. The economic inventory manage- 
ment model was used during this period. 

The model adopted in the 1970's has been revised systematically several times since 
its adoption by the Korean Air Force. However, the usefulness of this inventory model 
has decreased as the weapon systems used in the Korean Air Force have become more 
sophisticated, and the overall size of the Korean Air Force has increased. The annual 
budget for repairables takes approximately 60-70 percent of the total annual stock fund 
budget of the Korean Air Force, however the essentiality of high quality repairables in- 
ventory management was not realized. [Ref. 1: p. 1] 

In this chapter, the Korean Air Force Logistics system will be described. The 
chapter specifically deals with repairables management, organization, system parame- 
ters, and mathematical models related to repairable item management. This chapter will 
provide the basis for the analysis and suggested improvements to the Korean Air Force 


repairables inventory management system which will be described later. 


B. THE KOREAN AIR FORCE LOGISTICS COMMAND ORGANIZATION 

The Korean Air Force Logistics Command (KAFLC) 15 responsible for managing 
the allocation of the logistics resources to support maximum Weapon system operational 
availability. It functions as an intermediate echelon command among the tactical units 
of the Korean Air Force under the policies and planning of Air Force Headquarters. 
Figure 1 shows the organization of the KAFLC. The following paragraphs will explain 
the basic missions and responsibilities of the six major divisions of the command. 

The DMM (Directorate of Materials Management) is the focal point of material 
management for the Korean Air Force. It procures all materials according to its esti- 
mates of requirements, distributes the material to all of the tactical units and the sup- 
porting units under the policies of the KAFLC. Thus, the DMM is the equivalent of a 
wholesale level Inventory Control Point (ICP). The DMM manages approximately 


210,000 items which are allocated among the hundreds of item managers. To manage 
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rigure 1. The Organization of the Korean Air Force Logistics Command 


emcee), DM relics heavily on the Electronic Data Processing Center's (LDPC) 
computer services to access historical transaction data, compute demand forecasts, cur- 
rent asset data, etc. 

The DST (Depot of Storage and Transportation) is the centralized warehouse for 
the Korean Air Force where all procured and repaired materials, including consuniables, 
are stored. They control the physical inventory of the items needed by the ultunate user. 
The DST also deals with the disposal of sulwage, scrap, excess and obsolete material. 
The other important role of the DST is transporting all materials to arrive at the tactical 
unit when required. 

The DME (Depot of Maintenance and Equipment) iS thc in-house, depot level 
maintenance organization for the Korean Air Force. The primary task of the DME is 
the overhaul of aircraft on a scheduled basis (preventive maintenance). 1t performs the 
maintenance of aviation end items which are beyond the capabilities of the the interme- 
diate and organizational maintenance level facilitics. These items woukl be the major 
end items of an aircralt such as the engine, fuselage, gearbox and the direct ground 


support equipment for the operation of the aircraft. 


The DMEC (Depot of Maintenance, Electronics and Communication) performs the 
maintenance of radar equipment and its subassemblies, ground communication equip- 
ment and weather forecasting equipment. The DMEC also performs preventive main- 
tenance. 

The DMA (Depot of Maintenance and Ammunition) performs the maintenance of 
the precision measurement equipment (PME), airborne equipment and armament of the 
aircraft, such as laser guidance system. 

The EDPC (Electronic Data Processing Center) provides logistics software devel- 
opment and support, establishes job processsing standardization, collection of data re- 
lated to preventive and corrective maintenance, provides analysis and guidance for each 
of the other major divisions. The major objective of the EDPC is to assist each division 
in performing their specialized logistics missions within the Korean Air Force Logistics 


Command. 


C. REPAIRABLE ITEMS MANAGEMENT SYSTEM 
1. The General System Overview 

In the Korean Air Force an item of supply is designated as a repairable if it can 
be repaired faster and less expensively than it can be procured. Weapon systems in- 
stalled in aircraft have become increasingly sophisticated and complex. Hence, many 
weapon systems are made up of a number of subsystems which in turn are comprised 
of several repairable modules. Often the complexity of these individual modules are such 
that the personnel and equipment are not available at the end user level to repair failed 
units. Consequently, these modules are designated as repairables and failed units are 
returned to designated repair activities for repair. Therefore, repairables management 
has become an essential part of the Korean Air Force supply system in terms of decision 
making. Figure 2 shows the current Korean Air Force repairables management cycle. 

The Korean Air Force has two kinds of maintenance operations for repairables, 
the base maintenance (organizational and intermediate level) and the depot level main- 
tenance (DME, DMEC, DMA). It should be noted that not all depot level repairables 
are repaired at the Korean Air Force facilities. Because of the complexity, some repairs 
are done by the USAF under Foreign Military Sales (FMS) using commercial contrac- 
tors. As depicted in Figure 2, three kinds of repairables flows are shown. The first is the 
flow of carcasses to base level maintenance. The second is the flow of carcasses to depot 
level maintenance. The third level is the flow of carcasses from the DST to the USAF 
(the activities that coordinate with USAFLC for repair and procurement follow FMS 
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procedures). Note that the third flow implies the exchange of a carcass for a new item 
irom the SATE 

When a failure of a svstem or subsvstem occurs, the equipment is inspected to 
determine the cause of the failure. After isolation of the failure to a component or major 
assembly, further inspection is done to determine if base level repair is possible. In case 
of the failure of a major end item, an RFI (Ready-for-Issue) item from the base supply 
1$ used to replace the failed item and the aircraft returns to operational status. This in- 
sures the operational availability of an aircraft. If possible, the failed item 1s repaired 
by the Base Maintenance Squadron (organization or intermediate level of repair). After 
such a repair, the item is released to the Base Supply Squadron in serviceable condition. 

For the determination of the maintenance level, the Korean Air Force has a 
special code assigned to each replaceable component of a major end item. Thus, the 
level of maintenance for each component is predetermined by this code. The code is as 
follows: Ве. 


First code position 


X expendable item 
N repairable item 

Second code position 
D depot level maintenance item 
F field level maintenance item (intermediate) 
B base level maintenance item 

Third code position 
1 depot level maintenance 
D intermediate level maintenance 


3 no maintenance action 


After the Base Supply Squadron has issued a serviceable item from its stocks 
and when base level maintenance has been determined to be impossible, the Base Supply 
Squadron sends the failed item to the DST. Of course, all these actions are reported to 
and directed centrally by the DMM. The failed items from each base are placed in the 
DST facility and await repair. These carcasses are turned over to the depot maintenance 
facility (DME, DMEC, DMA) in batches under the approval of the DMM. In other 


words, all carcasses from each base are stocked in the DST facility during the quarter. 
At the end of each quarter. all carcasses are sent to the depot maintenance facility 
(DME. DMEC. DMA). However, at this time, item managers of the DMM consider 
the capacity limit of the depot maintenance facility.! 
2. Operating Characteristics of Repairables Management 

As mentioned earlier, there are several depot maintenance organizations in the 
Korean Air Force. However. their concept of operation and relations with the DMM 
are identical for the DMfE, DMEC and DMA. So, only the operation of DME will be 
covered here. The maintenance organization of the Korean Air Force is restricted in 
terms of its capacity to process incoming repairs and in terms of the level of technology 
to deal with the repair of complex systems. Thus, the MRS (Material Repair Schedule) 
and MRRL (Material Repair Return List) are established to manage these restrictions. 

a. Material Repair Schedule ( MRS) 

The MRS 15 the repair schedule for the depot maintenance facility. [t de- 
scribes the tvpe and quantity of repairable material which can be repaired during the 
next fiscal year by repair depots operated by the Korean Air Force. In September of the 
fiscal year, all item managers estimate the demand for each repairable item for the next 
year, and the total repair quantity, based on data from the last several vears. The inte- 
grated results for each item are provided to the DME, DMEC, and DMA. Then, the 
DME, DMEC and DMA set up their repair schedules after discussion and coordination 
with DMM. The forecasted maintenance requirements which are not covered by the 
MRS are turned over to the MRRL. The repair quantity for the MRS 15 set up on a 
quarterly basis and carcasses are issued from the DST based on capacity of the DME. 

b. Material Repair Return List (MRRL) 

This is the list of repairables for which the Korean Air Force does not have 
depot level maintenance capability, either due to technology or capacity restrictions. In 
the case of an MRRL repair, the DST sends the carcasses to a USAF depot facility. 
Upon the receipt of the carcass from the Korean Air Force, the USAF returns a serv- 
iceable unit to the DST. Serviceable units from both the USAF and Korean Air Force 
depot facility are integrated at the DST to make up the serviceable stocks which are 


available to support for the tactical unit. 


1 This capacity limit established based on the MRS, see paragraph 2.a. above. 


c. Ihe Measure of Performance 

The Korean Air Force supply system (wholesale level) uses two measures 
of performance. The first one is fill rate. It is also used by the base level supply system. 
The second measure is supplv response time. It measures the length of time elasped for 
base backorders to be satisfied bv the stock from the DST. In other words, it 1s the 
length of time from the placement of an order with the DMM to the receipt of the or- 
dered item. Notice that the DMM uses this more as a priority rule for issuing the stock 
than as a general measure of performance for management of the repairable items. 
Currently, the Korean Air Force assumes that any backorders for spares result in an 
aircraft being “not operationally ready because of supply” (NORS). Thus, the Korean 
Air Force has developed codes which indicate the NORS condition. Each code provides 
the maximum supply response time requirement to fill backorders. Table 1 summarizes 
the supply response time requirements for each demand with priority and. Table 2 de- 


scribes the meaning of each NORS code. [Ref. 4: p. 81] 


Table 1. SUPPLY RESPONSE TIME REQUIREMENTS 


Tepe of order Response time (day) 
G, K NORS 


A, F NORS 





Table 2. DESCRIPTION OF NORS CODE 
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D. MATHEMATICAL BACKGROUND OF REPAIRABLES MANAGEMENT 
1. DDR (Daily Demand Rate) 
This is the average number of demands per day and it is computed following 
formula: 


365 
DDR- 2 
-— 365 

i=1 


where: 


D, = demand on the ith dav of the year 


2. DRP (Depot Repair Percent) 
It 1s computed in the following manner for each repairable item from the last 


three vears of historical data: 


ETS 


И СОО 


х 100 


where: 
RTS - repair this station 
NRTS = non-repair this station 
COND = condemned 


RTS refers to the number of units repaired at the Korean Air Force depot. 
Conversely NRTS refers to the number of units repaired by the commercial contractors 
or by the USAF under FMS policies. Thus, NDRP (Non-Depot Repair Percentage) 1s 
the complement of DRP. 

3. RCT (Repair Cycle Time) 

The Korean Air Force has two kinds of repair cycle time. The one is the repair 
cycle time associated with base maintenance and the other is associated with the depot 
level maintenance. Usually RCT refers to the time allowed for the depot level mainte- 
nance only. This time allowance is the policy of the Korean Air Force. So, RCT is 


constrained3 to be not less than 30 days and not more than 120 days. 


2 This is computed by the end of year. The time period of this computation starts from the 
first day of the year to the end of the year. 


3 This is the policy of the Korean Air Force, not a value which is computed from actual 
observations which are recorded in the inventory management database. 


4. RCQ (Repair Cycle Quantity) 
This refers to the quantitv of the repairable item which must be stocked to meet 
demands during the repair cycle time. RCQ 1s applied to MRS items only and computed 


bv the following formula: 


КСО = DDK x ФЕР ХАКЕР 


where: 
DDR = daily demand rate 
DRP = depot repair percentage 
RCT = repair cycle time (days) 


5. OLQ (Operational Level Quantity) 
This is the stockage quantity which is required to cover forecasted demand 
during the time an MRRL item is in the hands of the USAF for repair. OLQ 1s related 
to MRRL as RCQ is to an MRS item. The formula for OLQ 1s as follow: 


OLQ = DDR x 60 day 


where: 


DDR = daily demand rate 


6. OST (Order and Shipping Time) 

OST is defined as the time elasped from the initiation of a procurement or re- 
pair order, to the receipt of the order. In the case of a procurement, the Korean Air 
Force constrains the OST to be not less than 120 days and not greater than 365 days. 
For MRRL repairables, OST 1s constrained to be not less than 220 days and not greater 
than 465 davs. In both cases, the upper and lower intervals are adopted to avoid an 
inventory stockout. The 100 days increment in the OST is due to the additional trans- 
portation time for MRRL items from Korea to the continental U.S.A. 

7. OSTQ (Order and Shipping Time Quantity) 
OSTQ refers to the quantity of repairable items needed to meet the demand 


during the order and shipping time. The formula for OSTQ is as follows: 


for MRS items, OSTO = DDR x NDRP x OST 
for MRRL items, 0510 - DDR x OST 


where: 


DDR = daily demand rate 
NDRP - non-depot repair percentage 


OST = order and shipping time 


8. SLQ (Safety Level Quantity) 
This refers to the quantity of repairables stocked to meet demand during delays 


in shipping, delays in maintenance, or unexpected increases in demand. The formula for 
SLQ is as follows: 


for MRS items, SLO = 3 х (АСО + OSTO) 
Ip MERE gems,  SLO = ,/3 х 0570 


Where: 
RCQ = repair cycle quantity 


OSTQ = order and shipping time quantity 


9. RO (Requisition Objective) 

Even with the assistance of the existing computer system, it is hard for the item 
managers to manage the more than 2,000 items (per item manager) that are assigned to 
them. Thus, the DMM designates items that have two or more requisitions per year as 
requisition objective items. The item managers will monitor these items more closely. 

Currently, the DMM applies a periodic inventory review model which 1s based 
on a policy of reviewing and repairing at fixed regular intervals to bring inventory levels 
back up to the Requisition Objective (RO). Repair orders are placed at predetermined 
intervals+ (the expected demand between intervals plus some allowance for the variabil- 
itv of demand must be considered). 

The requisitioning objective consists of three terms. These are the order and 
shipping time quantitv (OSTQ), the repair cycle quantity (RCQ) and safety level quan- 
tity (SLQ). The requisitioning objective is the sum of these components. In Figure 3, 
the repair order quantity is simply the difference between RO and the on hand inventory 
at time of review. 

Since the system parameters such as OST, RCT, DDR, DRP and NDRP differ 


at each review period, the level of the RO may also be different at the time of each re- 


4 This predetermined interval is the repair induction interval on the MRS. The current length 
of the interval 1s 2 months. 
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Fieure 3. The Korean Air Force Repairables Inventory Periodic Review 


view. Inventory reviews are usually done every 2 montis. Ulie formula jor XO 15 БЛОХ 
below 


for MRS items, RO = SLO + RCO + 05ТО 


where: 
SLQ = 4/3 x (hee Ose 
RCQ = DDR x DRP x RCT 


OSIQ = DDR x NDRP x OST 


Гог УЛ 1 TIEN KO = ООО о О 


Where: 
SLQ = 3 XOS? 
OLQ = DDR х 60 days 


OSTQ = DDR x OST 


10. History of RO Computation 
This section shows the history of the RO computation formula. The Korean 
Air Force tries to determine optimal stock level for the repairables inventory manage- 
ment svstem. However, the results of changing the RO formulation without adequate 
consideration of the svstem parameters may result in excess stock for some items or 
stockouts for some items. 
a. Before 1973 


RO - OLQ + SLO + OSTO 


where: 
OLQ = DDR x 75 
SLQ = DDR x 15 
OSTQ = DDR x 30 


b. Between 1973 and 1976 


For MRS, RO = DDR x 120 
For MRRL, RO = RKO + SLO +F OSTO 


where: 
TE DDR x 15 
80-13 х (RCQ + OSTO) 
OSTQ = DDR x $4 


c. Between 1976 and 1979 


RO We уо жо шо 


whare: 
NGO = DDR x DRP=x RET 
SEQ —9RCQ --O0STQ t NCQ) 
NGO = DDR x NDRP x NCT 
@s70 = DDR x NDRP x OST 
NCT = NRTS / Condemn Times 


5 This (NCT) is a time spent for decisioning a condemnation among NRTS items. 


d. Between 1979 and 1984 
RO = RCQ + SLQ + OSIQ 


where: 
RCO) = DDR rie 
5,О- СО 0 Р 7 р тоот га. 
OSTQ = DDR x NDRP x ROST 


e. After 1984 


For MRS, RO = RCQ + SLQ + OSTQ 


where: 
RCO = DDR ОКЕ ЕЕ 
SLQ = /3 x (RCQ + OSTQ) 
OSTO = DDR >R DREE TOT 


For MRRL, ВО - OLQ 4 SLQ 4 OSTQ 


Maere. 
ОГО = DDR x 60 
SLQ = /3 OSTQ 
OSTQ = DDR x OST 


E. SUMMARY AND PROBLEM STATEMENT 

As mentioned earlier, military inventory managers tend to want to increase their 
operational availability rather than minimize inventory costs. However, in the case of 
the Korean Air Force repairables management process, the mathematical model is in- 
sufficient to adequately model the stochastic characteristics of the inventory manage- 
ment problem. In the RO (Requisition Objective) computation it does not specify the 
wear-out rate or the regeneration rate, which are important parameters in dealing with 
repairables inventory management. 

The DDR (Daily Demand Rate) computation is the mean of a stochastic process. 
However, the current inventory model treats demand as if it were deterministic. In the 
formula for SLQ (Safety Stock Quantity), they considered the stochastic nature of de- 


mand but the actual probability distribution is not considered at all. As evidence of the 
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models shortcomings, the Korean Air Force repairables management tends to stock 


large quantiues of spare parts and components in order to avoid stockout situations. 
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HI. OVERVIEW OF REPAIRABLE MANAGEMENT MODEL 


A. INTRODUCTION 

Inventorv svstems such as the one described in the previous chapter can be classified 
rather broadly as a multiechelon system with repair. Analytical studies of multiechelon 
svstems have shown the computation in such models to be hard, so either simplifving 
assumptions or approximations are necessary. The introduction of the repair aspect into 
the svstem certainlv complicates the structure of a model even at single echelon levels. 

Suppose a repairable system, consisting of one Inventory Control Point (ICP), one 
stock point and one overhaul and repair activitv controls the inventory of a single item. 
Demands from various tactical units are placed only at the stock point. The system has 
continuous updating of records, 1.e., transaction reporting. When items wear out or fail, 
the tactical unit can either scrap the item or return it to the depot repair facility (DME, 
DMA, DMEC). After inspection, the depot can either scrap the item or repair and re- 
turn it to the RFI inventory. 

In this chapter, we will discuss various repairable inventory models based on two 
distinguishing characteristics. One is a deterministic model, the other is a stochastic 
model. In case of deterministic model, we will describe two models, “continuous sup- 
plement” model and “substitution” model. For the stochastic case, we will describe five 
models based on current usage in U.S military. Each description will include model 


formulation and solution techniques. 


B. DETERMINISTIC INVENTORY MODEL 

These deterministic inventory models were developed by David A. Schrady and 
others between 1966 and 1967. These deterministic inventory models are described in a 
research report for the U.S Navy repairable inventory management system [Ref. 5, 6, 
7]. The first model calls for the repair activity to induct a batch of carcasses when the 
NRFI inventory level reaches a certain point. With a deterministic system, this rule in- 
sures regularly spaced repair inductions of fixed batch size, which is a simple process for 
the repair activity to manage. Repair batch sizes can be maintained if there 1s contin- 
uous supplementing of repaired RFI items with new procurement. The procurement 
order will be issued with enough lead time so that when the on-hand RFI stock drops 
to zero, a procurement quantity will arrive. The two inventories, RFI and NRFI, will 


have time histories as shown in Figure 4. Note that in this, the “continuous supplement” 
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model, the procurement trigger 1s in the RFI inventory and the repair trigger is in the 
NRFI inventory. [Ref. 6, 7: pp. 9-10, pp. 392-395] 

The second model is suggested by noticing that there is a trade-off between stock 
held in RFI condition and stock held in NRFI condition. NRFI carcasses are a legiti- 
mate resource because they require only repair to be restored to full usefulness. But the 
value of this resource is less than the value of the RFI resource by at least the cost of 
labor and replacement parts. Thus, if inventory is to be held in the system, it would be 
slightly cheaper to hold it in NRFI condition than in RFI condition. [Ref. 6, 7: p. 10, 
р. №5] 

Notice that the “substitution” policy minimizes the RFI inventory while the "con- 
tinuous supplement” policy minimizes the NRFI inventory [Ref. 6: p. 13]. The next 
section describes the mathematics of these two models. 

l. Analysis of The Continuous Supplement Policy Model 

а. Notation 
Under the assumption of deterministic demand and lead times, it is not 
necessary tO maintain a safety stock. Whenever the RFI inventory reaches zero, a pro- 
curement quantity will arrive. Thus, the system is never out of stock. Between pro- 
curement arrivals, depleting RFI stock is replenished with repaired items. Also, a 
repetitive system will be established regardless of the initial provisioning policy, so that 
it is sufficient to analvze only one cycle to determine system characteristics. Figure 4, 
shows basic characteristics of the “continuous supplement” policy model. The basic 
notation as follows [Ref. 5: pp. 12-22]: 
D = annual demand rate 
Г, 


Г, = repair cycle 


procurement cycle 


t — procurement lead time 


t, — repair lead time 


r, = field recovery rate 
r, ^ overhaul and repair recovery rate 


procured lot size 


t° 
II 


Q, = repaired lot size 
h, = holding cost for RFI 
h, = holding cost for NRFI 


reorder point for procurement 


эс 
l 


A, = reorder point for repair 

n = number of repair cycles per procurement cycle, и = Г/Т, 
A, = fixed procurement cost per one cycle 

4, 


fixed repair cost per induction 


b. Model Formulation 

The objective of the model is to determine an optimal procurement quan- 
tity Q; and an optimal repair quantity Q; that minimize the average annual cost. These 
optimal values, coupled with respective reorder points, X, and X,, constitute an optimal 
operating doctrine. To determine the average annual variable cost, the costs per repair 
reorder cycle must be investigated. The product of per cycle costs and the number of 
cycles per vear will then vield the average annual cost. The total cost for any given cycle 
T, 1s the sum of procurement, repair and holding costs. 

Since there 1s only one procurement per reorder cycle, the variable order 
cost is the actual cost of the items ordered and can be expressed as Q,C,. The fixed 
procurement cost for one cycle is A, as mentioned. To determine the repair costs per 
cycle it is necessary to compute the number of repair cycles per procurement cycle ( 


п = T,/7T;, nis an integer). In Figure 4: 


1, = I + Г, 
where: 
n = KT, , O<K<1 
DO U 


The cost of items repaired per procurement cycle is C,Q, and the fixed repair cost 1s 
An, Where A, 1s defined to be the fixed repair cost per induction. The holding costs per 
procurement cycle will be h,A;, + A Ar, where Az is the area under RFI curve and 
Ar, is the area under the NRFI curve. To compute the area under RFI curve, consider 
Figure 5 showing RFI inventory for one procurement cycle. Since D is constant and 


known, the area of U, Aj, is given by: 


I 
Ау - > A, o) 
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Figure 4. Repair Cycle and System for “Continuous Supplement” Model 


where: 
1, = K ic 
C= 01 


The formula will reduce to: 
Ау = (10 -— D i) 2 


Since n is an integer, the area of V is the sum of n-1 trapezoids each having 


one side of length 7,. To determine the area of each trapezoid, a recursive relation 1s 
developed. let: 


Ay 


! 


area of the ith trapezoid in V (see Figure 4) 


== Á y, + Ay, + ..... ти 


where: 


1 
Аи =5 1420, % 2(Q@, — a D) — DT] 


Ay, = чу 1140 + 2(0,- ,0)- 3074 


therefore, the general formula for the area of V: 


Ay, = 5 14210, 4 2(Q - 4D) - Qi- DT]. i= 12... a= 1. 


n-i 
So the area of V is, 4, — Y Ay: 
i=] 


n—1 n—1 


Ay = ч ТПО о Пе | 


(El i=l 


x TQ: n(n = 1) + 2An=1Q, = An= )ар — DT wn 1) + ртн – 1) 


1 
5 Ти - 1, - &5) + (т - n)Q, — (n — 1)0Т) 
Finally, area of W in Figure 4, Ay 1s easily calculated as follow: 


| 
Ay = pu 
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where: 
А — ЮЛ, 
B O DBDs DBE 


Thus, the area of W: 
] 
Ay = 5 iuo meer -— 0 DAPI 
The total area under RFI curve is: 


Аг = Ay + Ay + Ay 


= nO, = — Dá + + T,(2(n— (Q) — Dn + ( — n)Q, — (т — mDT;) 
Since n is an integer and the build up rate of NRFI items, „0, 15 constant, 
the area under the NRFI curve is simply n times the area under the repair cycle curve 
БҮЙ 2». Sothe NKFI holding cost per procurement cvcle is 27,0; / 2r. The total 
cost per procurement cycle becomes: 








š A421 C2021 7,0, Т, 
Кт -06С-444 ur Яг Бонн + mcum £ EU 
The total average annual cost is then: 
12 Кт, - ОС A К СО; hQ, 4 һ Ат, 
Т, Т, Т, T, L 2r. Т 
Simplifving the total average annual cost: 
A,RD h h h 
K = + C, RD + А,” + 22 яг ле = АКО, ZI по 
Qi 2r, 2 2 


where: 
1, - 0,170:9 
ВЕ, n- L2,5,.. 
ЕКО р ког, Оче | 
1 + L= I, 
ЕО yD = Q, | RD 
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R = ] ын Vola 


Note that the terms CRD and CoD are independento Omand O asme neo 
affect the operating doctrine. Therefore it 1s appropriate to redefine the average annual 
cost of ordering, repairing and holding as follow: 








A RD Aror D hQ, һ 0, h Q 
О, F О; F 2r ES 2 — 11190 = 5 


From differential calculus the optimal Q, and Q, will be those that satisfy the equations, 











СК eK 

= = бе 0 

СО 00, 

Taking the partial derivatives results in: 

ак A RD Л, 

ПО ИС 2 
AK А”) h h 
a = eee — kh, + — 
СО, QE 2r, 2 


Solving the equation to get the optimal] procurement and repair quantities, 


: ARD 2А,гу 0 
Q, = h. d СЕО) 
1 2 238] 


2. The Substitution Policy Model 
a. Notation 
Notice that this model minimizes the RFI inventorv. Figure 5 shows basic 

characteristics of the "substitution " policy model. The basic notation as follows [Ref. 
6, 7: pp. 13-17, pp. 393-396]: 

О» = procurement quantity 

Оһ = repair batch size 

d — demand rate 

(1 —r) = scrap rate, r is the recovery rate measured 

tp, tp = procurement and repair lead times 


Ap = fixed procurement cost per order 
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A, = fixed repair batch induction cost per batch 
RFI holding cost 
h, = NRFI holding cost 


e 
ee 


И = system cycle time 


b. Model Formulation 
Referring to Figure 5, define T, as the time period during which inductions 
are suspended and the repair is simply accumulating NRFI carcasses. It may be deter- 


mined that: 
Та = (Qp + Ор) 1а 


Next. let n be the number of inductions per cvcle. In general, it will not be possible to 
insure that the last induction before T, begins will reduce the NRFI stock to zero. The 
residual NRFI stock when 7, begins will be some fraction of the net loss in NRFI items 
per repair cycle, where the repair cycle is the time between regular consecutive in- 
ducuons. The net loss of NRFI inventory over the period between successive in- 


ductions: 
шэг S ОБЕ) 


Thus, the residual when T, begins will be some fraction of Q,(1 —r), call it BQ,(1 — 7) 
Where 0< = 1. While the f factor is unavoidable due to the requirement that there be 
an integral number of repair batches in each cycle, we shall set 6 = 0 in subsequent de- 
velopments. 

Now the first batch after inductions are resumed takes the amount O, from 
the NRFI inventory. Subsequent inductions cause a net reduction of only О„( -- r) 
items. Thus, the amount of NRFI items available for the (n-1) inductions are 
ат, — QO, or On{(r —1) + rQ,. Dividing the last expression by Q,(1 — r) yield: 


= ИЕ (———)00% 
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Figure 5. Repair Cycles and System for “Substitution” Model 
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The system cycle ume is given by: 


Qp 


— 


T= пов + Ч _ 
Ë d ET 


where: 


n = ( - )0:Qs. 


| си 





From the Figure 7, the area under the RFI curve, over one complete cycle A,, is then: 


Ay = )سے‎ 0.9, + (“FOF 





Referring again to Figure 7 the area the NRFI curve: 


] 
А; = (42-2100 + ОР] 


+ 


Now we get total cost per unit times expression: 








TC a(1-r) 4 prd hor ee 
7 = А0 + Q + 777 008 + (7—00 + 10 + Or 
where: 
DM aa 
d di =— م‎ 


The total cost under substitution model: 
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The optimal order quantities are obtained by setting the partial deriatives, 


mae respect tO GO, and ОЈ: 


TC т 
Qp COR 
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Next we get new formula for optimal order quantities: 
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The optimal order quantities for repair and procurement: 


. 2A4pd(l-— r) * 24„4 
Op = a ee 202. 5-1 
ИШ о h, + h, 


3. Summary and Evaluation 

An important result of deterministic inventory model analysis is the realization 
that economic order quantities do not vary directly with demand but as the square root 
of demand. This nonlinear relationship could explain why many organizations experi- 
ence inventory problems. Intuitive inventory policies are undermined by difficult-to- 
grasp nonlinearities. [Ref. 8: p. 141] 

Although the models described in this section are evidently applicable to re- 
pairables management, thev are fraught with many limitations. The reasonableness of 
model assumptions and their sensitivity to actual conditions determines the utility of anv 
particular model. For example, the two models assume that the demand for a repairable 
item is known with certainty. Also stockouts either are not permitted to exist or are 
backordered and satisfied when replenishments are received. 

In terms of the Korean Air Force repairable management system, the models 
described here are difficult to adapt to the Korean Air Force system. This is because 
of demand variability, relatively long procurement lead time, repair capacity and capa- 
bilitv, a number of backorder situations, and current system management concept. The 
other important point is that the Korean Air Force often encounters NORS conditions 
which results in decreased operational availability. However, these two deterministic 
models do not deal with stockouts. These reasons didn’t allow to use for the Korean 


Air Force repairable management system. 


C. STOCHASTIC INVENTORY MODELS 
As the weapon systems installed in modern aircraft become increasingly sophisti- 
cated and complex, repairable items represents an important subset of the total inven- 


tory of items which are managed by the military supply system. 
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In general, repairable items are supported bv a two-echelon inventory and repair 
system. When a repairable item fails at the base level, it is returned to base supply and 
a new serviceable unit is issued from extra RFI stock of the base. If possible, the failed 
item is then repaired by the base maintenance organization and then stored in RFI 
conditions at base supplv. Sometimes, however, the failed item must be returned to the 
repair depot where more sophisticated equipment and specialized skills are available to 
repair it. 

When we consider the condemnation of repairables, there should be an inflow of 
newlv manufactured items to depot level supply should look at the material flow of all 
echelons, and then, place a replenishment order to resupply the condemned repairables. 
At times, forward base locations will be supplied from another closely located base. This 
is called lateral resupply. For other items, a manufacturer may provide both the source 
of procurement for new assets and the source of repair for failed items. 

The decisions concerning which maintenance levels will repair the failed items sub- 
sequentlv affects the supplv support provided at the organizational level. If an item is 
repaired at the organizational level, repair equipment and maintenance personnel must 
be made available at that level. However, if the item 1s not repairable at that level, then 
the question is whether the item can be replaced at the organizational level. Organiza- 
tional level repairables will normally be transferred to the next higher echelon of repair 
if the repairs cannot be accomplished at the organizational level. Similarly, repairables 
which cannot be repaired by the intermediate level are usually sent to the depot level for 
repair. 

The difference between stochastic models and deterministic models is that if de- 
mands or future requirements are uncertain, then regardless of the stockage policy 
adopted there is generally some probability that available stock levels will be insufficient 
to meet demand [Ref. 9: pp. 187-188]. In this section, we will describe five stochastic 
models based on current usage. 

]. Model Analysis 

a. Analysis of METRIC Model 
(1, Background. The METRIC model was developed over a period of 
vears bv a research group at the RAND corporation and presented in the literature by 
Sherbrooke (1968) and extened by Muckstadt (1973). METRIC was developed with the 
ultimate goal of implementation and a slightly modified version of METRIC was actu- 
ally implemented by the USAF. [Ref. 10. 11, 12] 
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METRIC 15 а mathematical model of a base-depot supplv svstem in 
which item demand 1s assumed to be compound Poisson with a mean value estimated 
by a Bayesian procedure. The model is also capable of determining base and depot stock 
levels for stock levels of a group of recoverable items. Its governing purpose 1s to opti- 
mize system performance for specified levels of system investment. It is designed for 
application at the weapon-system level, where a particular line item may be demanded 
at several bases and the bases are supported by one central depot. The support depot 
may vary by item as in the item-manager system or it may be fixed as in the weapon- 
system storage site concept. The general purposes of this model are optimization, re- 
distribution and evaluation. [Ref. 13: p. 2] 

The basic METRIC model considers a two-echelon system in which 
independent bases (lower echelon) are supported bv a repair depot (upper echelon). At 
the occurrence of a failure (it could be more than one), the failed item is either replaced 
by available base stock or back ordered if the base stock 1s not available. The item 1s 
inspected to determine the extent of repair required. If the repair can be made at the 
base, the unrepaired item enters base repair. If the item cannot be repaired at the base 
level, it is shipped to the depot. Upon shipping the item to the depot, the base places 
an order with the depot for a replacement, so that the inventory position for item 1 at 
base } can be maintained. 

(2) Mathematical Assumptions 


1. System Objective of Minimizing the Expected Number of Backorders: The objective 
will be to minimize the sum of expected backorders on all recoverable items at all 
bases pertinent to a specific Weapon system. Thus, unless all bases are identical, 
the expected number of backorders will vary by base. Expected backorders take a 
fixed period of time and add together the number of days on which any unit of any 
item at any base is backordered. In order to minimize the expected number of 
backorders one divides this number by the length of the period and taking the ex- 
pected value. This vields a number that is independent of period length. This 15 
the value we seek to minimize. [Ref. 14: pp. 126-131] 


2. Compound Poisson Demand: We assume that demand for each item 1s described by 
a logarithmic Poisson process. Compound Poisson processes are generalizations 
of Poisson processes; the compound processes allow the flexibility of incorporating 
more parameters, yet retain the simple analvtical properties of the Poisson. The 
logarithmtic Poisson is obtained by considering batches of demands where the 
number of batches follows a Poisson process and the number of demands per batch 
has a logarithmic distribution. 


3. Demand is Stationary over the Prediction Period: It is assumed that the distribution 
of demand over some future period of interest 1s stationary. 
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4. Repair Decision Deponds on the Complexity: The assumption is that the decision to 
repair a unit at base level or send it to the depot is a function onlv of the type of 
malfunction and the base maintenance capability. 


сл 


Lateral Resupply is Ignored: When a unit is shipped from base to depot for repair, 
a serviceable replacement will be resupplied from the depot if possible. If the depot 
has no unit on the shelf, the base must wait until a unit emerges from depot repair. 


6. System is Conservative: Consider a particular stock item demanded at base j. We 
assume that a unit of stock has a probability r, of being repairable at base j, with 
a repair time drawn at random from a base repair distribution with mean 4; a 
probabihty (1 — rj) of being depot repairable, with an order and shipping-time 
distribution having mean O, and a depot repair distribution having mean D. This 
implies that there are no condemnations or that the svstem is conservative, as the 
name “recoverable item” suggests. A higher condemnation rate usually indicates 
that the item should be redesigned. The condemnation rate must be considered for 


procurement purposes, but the procurement process is not considered in the 
METRIC optimization. 


7. The Depot Does Not Batch Units of Recoverable Item for Repair: The model as- 
sumes that depot repair begins when the repairable base turn-in arrives at the de- 
pot. This appears to be a reasonable approximation to current depot scheduling 
practice. Since METRIC economizes by buying fewer high-cost items, these are 
the items which are most likely to be in short supply. 


$. Demand Data from Different Bases Can Be Pooled: We assume that demand from 
the several bases can be pooled in some manner so that a composite initial estimate 
of demand per flving hour can be obtained. The pooled estimate can be obtained 
by a simple averaging technique or a more sophisticated procedure such as expo- 
nentıal smoothing. METRIC multiples this number by the flying hours per month 
for each item. 

(3, Model Formulation. By mathematical assumptions, demand in this 
system follows a compound Poisson process. A compound Poisson process may be 
thought of as a series of customers who arrive following a Poisson process. each of 
whom can demand an amount that is independently distributed according to a com- 
pounding distribution. Assume that item 1 is stocked at each of base j, and the cus- 
tomers who place demand for the item at each base have a known mean arrival rate of 
/,, }= 1,2,3,......,J. When a customer arrives at base to place one or several demands, he 
turns in an equal number of carcasses. These carcasses can be repaired at base level with 
probability of r,, while (1 — r,) 15 the probability that they must be repaired at the de- 
pot. The arrival of carcasses from base j at the depot is described by a Poisson process 
whose mean is (1 — r,) times the mean of the Poisson customer arrival process at the 
base j. Therefore, the total demand at the depot for item i is compound Poisson, with 


mean customer arrival rate: 


3l 
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j=) 


where: 


Àj — mean arrival rate fot item 1 at base j 


Let /, be the mean demand per customer at base j. Then, the mean 


depot demand rate for unit i is: 


J J 
в = учат ба сп 
j=l j=l 

where: 


0, — the mean rate for item i at base j 


In the special case of the logarithmic Poisson process, the probability 
that x customer demands are in the repair/supply process 1s negative-binomial with pa- 
rameters q and k (Note: K=AT/ Ing where / 1s mean customer arrival rate and T is 


average resupply time). 


Ca 1)“ 
.. .. — m |  —— лизна 
P(x | AT) (K + x — I) Dia 


x * 012.499 4] ü 


where: 
q =the variance to mean ratio of P(x | A,T;) 
K 2 ATA d (q — 1) 
Г, = average resupply time for a demand for item 1 at base j. 
А, = base repair cycle time 
D = depot repair cycle time 
O,, = the order and shipping time 


6(S,) = average depot delay 


Then: 
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Since 1t takes an average of D time units for an arrival to complete 
the repair process, the probability distribution of the number of units in the depot repair 
cycle is compound Poisson with mean of AD. Hence, the expected number of units back 


padered at the depot is: 


BO(So|4D) = 2 (x — Sg) P(x12D) 


X> So 


As mentioned in the beginning, the objective of METRIC is to min- 
imize the sum of backorders for all item 1 and for all bases j within a budget constraint. 
Thus, the METRIC problem vill be represented as follows: 


на 
minimize У Sj) 


Б-1 j=1 


7 
subject to eus C 


i=1 j=l 


where: 


the decision variables, stock for item 1 at base j 
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C = the total amount of budget available 


the cost of item 1 


су 
| 


S, ^ the depot stock for item 1 


(4, Solution Technique. The METRIC problem is solved by using the 
generalized Lagrange Multiplier method suggested by Fox and Landi [Ref. 15: pp. 
258-261]. Let O be a Lagrange Multiplier associated with the budget constraint. [he 


Lagrange function 1s written: 


7 ж 
> > воду 2919) = оу 272) 


i=1 J=1 i=1 j=1 


The auxiliary problem attempts to minimize this equation. By trial 


and error, we trv to find the value of ® which satisfies a given constraint. Therefore, 
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we need to solve the above equation for several values of «b, and choose that value of 
Ф for which the required resources are closest to the budget limit. Fox and Landi sug- 
gest a binary search procedure. Their computational experience was that at most six 
bisections Were required to obtain budget allocations that were within one half of 1% 
of the original budget C. 

The objective function is separable in the items. Dropping the sub- 


script 1 in the original problem allows us to rewrite the equation as: 


1 
шіп ),80(5/14/1) - ФСбу- ФСбо 
i=] 
Since BO(S,|4,T7,) is discretely convex for a given S,,, the eptimum 


base level 1s obtained by simply finding the smallest non-negative integer satisfying: 
BO(S; + 1|4;T)) — BO(S,|7;7)) = OC Sip 


(5) Model Evaluation. The current Korean Air Force Logistics system's 
objective is maximizing operational availability subject to a budget constraint. This is 
sinular to the METRIC model. The objective of minimizing the expected number of 
backorders is related to the objectives of the Korean Air Force Logistics system. How- 
ever, two assumptions are distinctly violated. One is batch size repair policy is not used, 
the other is lateral resupply 1s ignored. Currently, these assumptions are not used in the 
Korean Air Force because of relatively long procurement leadtimes and repair capability 
limitations. It may be possible for the Korean Air Force computer resources to handle 
the data requirements for the METRIC model, but it is doubtful for actual usage. 
However, the METRIC model is obviously related to the Korean Air Force in terms of 
its basic concept and system objective. 

b. Mod-METRIC Model 

(1j Background. Mod-METRIC model was developed by Muckstadt to 
deal with. problems which the METRIC approach did not consider. Mod-METRIC 
considers the relations in parts hierarchy and tries to solve this multi-indenture level 
problem. Mod-METRIC was implemented by the U.S. Air Force as the method for 
computing repairable stock levels for the F-15 weapon system. [Ref. 16: p. 472] 

Most repairables contain subassemblies or other components which 
are also repairable. For example an aircraft engine, it may have modules for intake, 


combustion and exhaust. If an engine fails, it is replaced by a serviceable engine from 
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base stock. The failed component is then repaired either at the base or depot level de- 
pending upon the complexity of repair. A serviceable module from the base stock, if 
available, will replace the failed module, and the repaired engine is placed in base engine 
stock. 

(2) Mathematical Assumptions. The assumptions stated in METRIC are 
also applicable to Mod-METRIC except that Mod-METRIC assumes that the demand 
process 1s the simple Poisson process. [Ref. 16 pp. 475-476] 

(3, Model Formulation. In METRIC, the objective is to minimize ex- 
pected backorders for all items subject to an investment constraint. The Mod-METRIC 
objective, however, is to minimize the backorders for the end item, subject to an invest- 
ment constraint on the total dollars allocated to the end item and its components [Ref. 
16, 17: p. 475, p. 39]. This difference is caused by the hierarchical maintenance re- 
lationship between the module and the end item. As an example, an engine backorder 
indicates that an aircraft 1s missing an engine and is unavailable to perform its flying 
mission. Modules, on the other hand, are used only to repair engines. A backorder for 
a module only delays the repair of an engine. The impact of module backorders and 
engine backorders is clearly not the same. 

Let T; denote the average engine resupply time at a base. By the na- 
ture of repairables, 7; depends on several factors. When it 1s repaired at base level, J; 
would be the time it takes to flow through the base maintenance system. If it is repaired 
at the depot, 7, consists of the time to place the depot order for a serviceable part and 
to receive the part from the depot, assuming a serviceable asset 1s on hand at the depot. 
However, when there are no serviceable assets on hand at the depot, an additional delay 


1s included in the resupply time: 
Тр = пВ; + (1 — "4; + 9(5,р)) 


where: 


r, = the probability an engine will be repaired at base 1 


B, — the average resupply time, given an engine is repaired at base 1 
A, = the average order and ship time for an engine at base 1 

S, = the depot engine stock 

D = the average depot repair time 


óS,D = the delay days per demand 
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0S,D can be derived in the following manner. The expected number 
of engines back ordered at the depot 1s: 


BO(S|AD) = ) (X — S,)P(AIAD) 


X> S, 


where: 
I 


20 NE 


{= | 
; ^ the daily engine removal rate at base 1 


NY 
Il 
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In other words, this expression is the expected number of units on 
which delay is being incurred at a random point in time. Dividing this expression by the 


expected number of demands per dav yields a statistic which has the dimension of delay 
days per demand: 


oS DESEE 


— expected backorders/expected daily demand 


Further, B, can be divided into two components. That is, Z, 1s equal 
to the average remove and replace time, given that the necessary module is available, 


plus the expected delay due to the unavailability of the module which 1s required to гс- 
pair the engine. Then: 


B; = Ri + А; 
where: 


= 


1 


the average repair time at base 1 if modules are available 


A,= the average delay in base engine repair due to unavailability of a needed module 


Further, assume that the engine should be repaired by the failure of 


module j. Then, A, 1s the expected delay in engine base repair time due to a back order 
on module j at basei. Thus: 


Ау - 2 (Ay — Sg POS | AyTy) | Ay 
X; > 5; 


where: 
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^, 7 average number of daily removals of module j at base i 


| 


stock level of module j at base i 
I, = average resupply time for module j at base i 
Then: 
vw efe 


r, = the probability that a failure isolated to module j will be repaired at base level 


7 


В, 
ғ 


average base repair time for module j at base i 


average order and ship time for module j at base i 


^4 


D 


average depot repair time for module | 
Y (X - S)PQC16,D)/ 0D, 


= 


0 = 

| Х> 5, 

S, = the stock level of module j at the depot 
8 = 


V 


14: 


211 ЕЕ n) 


7 1 
Consequently, the expected delav in engine repair at base i due to 


module unavailability 15: 


J 
А = Одд ДА 
j=1 


Thus, expression 7, represents all the svstem components including 
the depot engine and modules stock level, and the base module stock level. By using this 
relationship. the engine and module stock levels can be derived. Finally, the math- 


ematical statement of the problem 16: 


1 
min 2, 2, (X, — SPP AT 


i=1 X> S; 


I J zi 
subject to 2 (GS, + 2 С)5) % 2 65, % С,5,< С 
т j2i 


id 


where: 


5 не Лехе! оГ base 1 
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C, — unit cost of an engine (end item) 
CE UNI COST Of OE 


C = dollar budget lımıt 


/4j Solution Technique. Unfortunately, the equation above is not sepa- 
rable, because I, 1s a complex function of the S;, Thus, Muckstadt recommended that 
it should be broken down into two parts; the component subproblem and the end item 
subproblem [Ref. 16: p. 477]. Even after the break-down of the problems, however, it 
requires the solution of many subproblems each of which corresponds to a particular 
division of the available budget between components and end items. The solution pro- 
cedure 1s as follows. 

First allocate the budget C into C, and C, for components and end 
items repectively. Allocate C, among components so as to minimize the expected end 
tems repair delays summed over all bases subject to the budget constraint C,. This 


problem is mathematically represented as: 


subject to ) (OS + > 65)<С 


This problem can be solved using the METRIC technique. Given the 
result from this step, compute the average resupply time T, for the components of each 
base. Then, allocate the remaining budget C, so as to minimize the expected end item 
base backorders. The METRIC budget allocation procedure mav again be used. The 
steps provide a set of proposed stock levels upon a given allocation of the budget among 
the end items and components. These steps then are repeated several times using new 
values for C, and C, to establish the best allocation. 

(5) Model Evaluation. Basically, the Mod-METRIC model has a similar 
Structure to METRIC. All assumptions used in Mod-METRIC model are applicable to 
METRIC except for the compound Poisson demand. Іп other words, this model is the 
same as METRIC for the Korean Air Force repairable management system. The only 
disadvantage is that the Mod-METRIC model was developed for new weapon systems 


such as the F-15 aircraft. However, the Korean Air Force has a lot of non-standard 
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Weapon systems. These non-standard weapon systems are causing problems such as 
long procurement leadtimes. 
c. Dyna-METRIC Model 

(13 Background. | Dyna-METRIC was developed by the RAND Corpo- 
ration to provide an analytic method for studying the transient behaviour of 
component-repair and inventory systems under time-dependent operational demands 
and logistics decisions like those that might be experienced in wartime [Ref. 14, 18]. 
Note that the past work regarding the repairable item stockage prior to Dyna- METRIC 
only dealt with a steady state inventory system with constant average demand and ser- 
vice rates. These steady state assumptions may provide a good approximation during 
peace time operations. But in wartime demands for components may suddenly increase 
very much relative to the previous peacetime operation and then may decrease gradually 
Or, in some case, drastically due to attrition in the system. [Ref. 14: pp. 1-6] 

A kev characteristic of Dyna-METRIC is its ability to deal with the 
dvnamic or transient demands placed on component repair and inventory support 
caused bv time variables in a scenario that includes sortie rates, mission changes, phased 
arrival of component repair resources, interruptions of transportation, and the like. all 
Oreemich would be experienced in wartime. It computes how given resource levels and 
process times would contribute to the Wartime capability. By exploiting the mathemat- 
ical structures of its underlving equations, Dyna-METRIC suggests the alternative cost 
effective repair or stockage resource purchases which would achieve a target aircraft 
availability goal throughout the wartime scenario. 

Dyna-METRIC considers a three echelon inventory repair system. 
Each base has an in-house repair facility which may have various test and repair capa- 
bilities. This base repair facility may be supported by several Centralized Intermediate 
Repair Facilities (CIRFs). Each operating base is capable of doing only limited types 
of maintenance, usually limited to simple removal and replacement operations at the 
flight line. It should be noticed that some of the bases have a CIRF while others do not 
have direct flows of parts to the depot. A depot is represented as existing outside of the 
model. It is seen from the model's point of view, as an infinite source of supply located 
some order and ship time away. 

/2; Mathematical Assumptions. The major assumptions’ of 


Dyna-METRIC which distinguish it most from the others are shown below. 


1. Demand for items are generated by a nonhomogeneous Poisson process with in- 
tensity function m(t), and mean value function r(t)= fms) ds. Тһе functions, 
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m(t) and r(t). will be defined later. The t denotes an arbitrary timc. Thus, the de- 
mand process 1s dynamic. 


2. The repair process is independent of the arrival process and has slack repair ca- 
pacity so that cach repairable item demanding repair immediately receives services 
with average service time based on the function F(s,t). F(s,t} wall be defined jater. 


3. No more than one subcomponent can fail or be demanded in the repair of cach 
assembly. 

(3) Time Dependent Pipeline Equations. Since the major objective of the 
system is to avoid the loss of aircraft mission capability due to a shortage of correctly 
functioning components on the aircraft, it is necessary to compute the number of com- 
ponents awaiting repair, being repaired, being on the way to and from another echelon 
of repair, and partially repaired but awaiting spare parts. Each state is a pipeline seg- 
ment is which characterized by a delay time that arriving components must spend in the 
pipeline before exiting the segment. The model expands cach component's expected 
pipeline size into a complete probability distribution for the number of components 
currently undergoing repair and on order, so the probability distribution for all compo- 
nents can be combined to estimate aircraft availabilitv and sorties. [Ref. 14: pp. 7-23] 

Under the assumptions that the probability distribution of repair 
time is independent of the failure process, the average number of components in the re- 
pair pipeline will be: 
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where: 
d = average daily failure rate 


7 = average repair time 


With the further assumption that demand has a Poisson probability distribution, the 


probability that there are K components in the pipeline at any point in time would be: 
NTT Ка 
P( K in pipeline ) = „це ““Ј К! 


However, in Dyna-METRIC the demand is a function of time so that: 
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d(t) 7 (failure per flying hour) 
x(flving hour sorties at time t) 
x(number of sorties day per aircraft at time t) 
x(quantitv of the component on the aircraft) 


x(percentage of aircraft with the component) 


Also, in place of a constant average repair time, T, the dvnamic model uses the proba- 


bilitv that a repair started at time s is not completed at time t. That is: 
F(t,s)= Probability\(component entering at s is still in repair at t) 


= Probabilitv(repair time > t-s when started at s) 


The average number of components in the pipeline is derived by combining those two 
functions. Consider only those components that arrived in an interval of time, As. сеп- 


tered at time s. Then: 


A/(t,s) = d(s) F(t,s)As 


where: 
A^(t.s) ^ expected number of components in the repair pipeline 
d(s) = daily failure rate at time s 
F(t,s) = Probability that a component is not out of repair by time t 


As = interval of time centered at s 


If we assume that the number of failures arriving in the interval As 1s independent of the 
number of failures arriving in similar intervals centered at other time other than s and 


F(t) 1s independent of the probability distribution generating the demand rate, then: 


A(t) = 2 d(s) Ба А 


SSI 


Further assume that As is very small, so that: 
t 
A(t) = fas Flr sidas 
0 
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With the additional assumption that the component failure probabil- 
itv distribution is Poisson, A(z) 1s the mean of a time varying (nonhomogeneous) Poisson 


process. That is, the probability of K components in repair at time t is: 
Р(К) = A(0*e "9 | К! 


where: 


A(t) = | dis) F(t,s) ds 
0 


(4) Time Dependent Component Performance Measure. The component 
measures typically computed by the Dyna-METRIC model are [Ref 14: pp. 24-29]: 
R(t) = readv rate at time t 
FR(t)= fill rate at time t 
EB(t)= expected back orders 
VBD(t)= variance of the backorders 
D'T(t)7 average cumulative demands bv time t 


S(t)= supply level at time t 
The ready rate is given by: 


S(t) 


К) = >, PKI ACM) 
K=0 
Since the definition of the fill rate 1s the probability that a component 
will be available when a demand 1s placed, it is therefore the probability that demands 
have left at least one component available, that is, the sum of the probabilities of de- 


mands less than the stock level. Expected backorders are given by: 


ЕВ) = 2, (K - s0)PLKI XQ) 


К>5(1) 


S 
= 2() — s()  » (st) — K)P(K | 2(0) 


K=0 


For K greater than s(t), there will be backorders of (K — s(r)). The 
probability ofany demand level, that is K, is P(K/ 4(1)), and the expected value of the 
backorders 1s merely the product of the various values the backorders can take on times 


the probability of a demand at that given value. The variance in backorders is given by: 


VBG) = 2, (K - sy P(KT40)) — (ЕВОУ 


K>S(t) 


(3) Time Dependent Optimal Determination of Spare Parts. The fact 
that pipelines have time-dependent probability distributions means that the optimal mix 
of spare components at one point in time may not be the optimal mix at another. Thus, 
the approach to take is to compute, for each time period of interest, the marginal in- 
crease in spare parts to achieve a given capability over those already input or determined 
for a previous time. In determining the supply level, the model attempts to provide 
enough spare parts to give the desired confidence at the lowest cost at each point in time 


of interest. Thus. the objective function is to minimize the total cost of spare parts. 


Ie 


S, ^ the spare parts level for component 1 


С, = the unit cost of component 1 
a = the desired confidence level 
K, = the non-mission capable rate not to be exceeded 


P(K,,S) = the probability that the non-mission capable rate is less than K, 
given a stock level S 


I = types of repairables required on each aircraft 


Then, the problem to solve is: 


I 
minimize 22: 
i=] 
swWujecuspessp( KAS) >, 97Г- 5; 


== “¿O 


where: 


S, = the input stock level or previous time optimization stock level for component 1 
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Assuming complete cannibalization, P(K,,S) equals: 


I 
Р(К,5) - | | (ОК) 
i=) 


where: 
S; 2 92,К, 


MOR) T T RAT) 


P (K) = the probability of exactlv K failures of component 1 


Ihe necessary condition for the performance constraint to be met is to have, 
P(Q,K,) =o, for each 1. Then marginal analysis is used to determine the best mix of 
additional components to achieve the desired goal. This process proceeds by investing 
in one additional component at a time which is selected by finding the component that 
gives the largest increase in the logarithms of the confidence level at the lowest cost. 


That is, we determine: 
Ain P (K,,S) | C, 


where: 
Ан Р (А = и 
S' = (SS ОМ) 


The component for which supply 1s increased one unit is the one whose index solves; 
тах А, /и P(K,.S) / C,. This process continues Until the given confidence levels achieved: 
At this point, the resulting value of S is the efficient solution of the base stockage 
problem. 

(6) Model Evaluation. The Dyna-METRIC model 1s relatively hard to 
adapt for the Korean Air Force. Especially since this model assumes that each repair- 
able item needing repair immediately receives service. This does not occur in the Korean 
Air Force. Currently, the Korean Air Force applies a batch repair policy for repairables, 
except in case of NORS condition aircraft. This is a serious violation of the assumption 
of Dvna-METRIC. So, the Dyna-METRIC model is difficult to use for the Korean Air 
Force repairable management system without using a continuous repair policy. 

d. Mathematical Models Used In The Navy 
(1) Background. The Uniform Inventory Control Program (UICP) was 


developed in August 1965 to provide a standard svstem to be used at all U.S Naval 
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Supply Systems Command (NAVSUP) ICPs. The Fleet Material Support Office 
(FMSO) under the direction of NAVSUP is responsible for the svstem design, ADP 
analysis. programming and documentation of UICP. The development of the UICP 
formulas for inventory levels follows the approach used by Hadley and Whitin. [Ref 
13, 19] 

(2, Mathematical Assumptions 


l. It is a continuous review system. Wholesale inventory levels requirements and as- 
sets are known by the Inventory Control Point (ICP) at all times. [Ref. 13, 19: pp. 
162-165, Chapter 3 Appedix A] 


2. It is a steady state environment. The key characteristics of the items managed by 
the ICP are assumed to be constant over the forecast period. These characteristics 
include the forecasted average values and variances of the random variables of the 
rate of customer demand. procurement leadtime, production leadtime, depot repair 
times, depot repair survival rate and the rate of carcass returns. 


(22 


. An order for procurement or repair is placed when the assets reach the reorder level 
or the repair level. 


4. Customer demands and carcass returns do not occur in more than one unit per 
transaction. 


Cn 


The unit procurement cost or repair cost of an item 1s independent of the magni- 
tude of order quantity or repair quantity. 


6. The cost of a backorder and the time-weighted cost of a backorder can be accu- 
rately quantified. 


/. The reorder level and repair level are always non-negative. 


8. The cost to hold one unit of stock in the inventory is proportional to the unit cost 
of the item. 


9. No interaction exists among families of items or individual nonfamily items or 
both. Each family or nonfanuly item’s inventory levels requirements are calculated 
independently of those of families or nonfamily items. 


10. The optimal inventory levels are determined by minimizing an average annual 
variable costs equation composed of order cost plus holding costs plus shortage 
costs. 

(3) UICP Depot Level Repairables (DLRs) Procurement Model. The 
repairables procurement model starts with a total variable cost (TVC) equation which 
is to be minimized. A notable difference from the model for consumables is the inclusion 
of receipt of RFI assets from a repair process in the DLR model. 

(a) Ordering Cost 


Let: 


A = administrative cost to order 
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D = forecasted quarterly recurring demand 
B = forecasted quarterly regenerations 
Q = economic order quantitv 


Then, 4(D — B)/Q is the expected number of procurements per vear. Thus, the annual 
ordering cost is given by 4(D — B)A/Q. 
(b) Holding Cost 


Let: 
= percent cost per dollar of inventory held annually 


= item cost (replacement) 


I 

C 

L — procurement leadtime 
T = repair turnaround time 
R 


— repair level 
F(x) = the function of probability distribution of leadtime 


demand 


The mean demand during a procurement leadtime 1s given by: 
DT + (D — B) x (L — 7) 


Thus, holding cost is given by: 


IC + R-(D- Bx L+Bx T + | (x — R)E(x) dx] 


x>R 


ket: 
À — shortage cost per requisition short 
E = military essentiality weight 


F = quarterly requisition rate 


Then: 


(AED 2210 
Г (x — R) F(x) dx 


x>R 
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It equals the cost of the expected number of backorders in a year. The TVC equation 
is symbolized byv: 


IVC = (D — B)A JO 
- IC[QI2 + R — {(D — B)xL + BxT) + f (x — R)Fx) dx] 
- GE[4(D — B)IQ] F| D Је — R)F(x) dx 


Then, by setting dIVC / dQ = 0: 


Q= |[5(0 – В)(4 + 2EF | D) | (x — R)F(x) dx | IC] 


x>R 


Also, bx settuing dT VC / dR = 0, we find: 


QICD 


от 


XR 


Because the expression и: — R)F(x)dx is the cumulative 
distribution for leadtime demand, this is the quantity defined as RISK. Note that since 
Q and R are related, the reorder quantity Q cannot be solved independently. Thus, 


LICP approximates Q by using a variation of the economic order quantity formula: 


8(D — B)A 
o> IC 


Then, QICD | (QICD + 4À¿EF(D — B)} can be computed for RISK determination. 

(4, UICP Depot Level Repairables Repair Model. The repair model also 
Starts with a total variable cost equation viewed largely independently of the procure- 
ment problem. Note that its time horizon is the depot level repair turn around time. 
The total variable cost equation for the repair model is the sum of the order cost, holding 
cost and backorder cost. 


4] 


J or: 


Q, = economic repair quantity 


A, = repair administrative order cost 
C, = repair price 
R, = repair level 


Ех) = probability distribution of demand during during repair turnaround time 


A, = repair shortage cost 


Then: 
IVC - 4min(D,A)A, | Q, 
+ 15 + R, — (D XT) + / (х — R,)F,(x)dx} 


T A,(Àmin(D,B)] Q3F |] D с (x — К.)Р,(х)ах 


By setting ol KC O aS 


О. = | 8 min (D,B)A, + ),EF| D | ЕЕЕ 
х> К, 


Adso, by settine ОЙ С 2017-1858 
| Е(х)ах = Q,IC,D | Q,IC,D + 41,EF mir(D,B) 
X20 


Again, Q, and R, are related. The UICP model approximates Q, as follows: 


О, = 6 min (р,В) 4, / IC, 


Then, Q,IC;D / Q,IC,D + 4¿,EFB can be computed by using the above results. 
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/5) Integrated Repairables Model. As stated earlier, the requirements 
computed by the procurement and repair models are accomplished independently of 
each other. As a result, this leads to a carcass constrained situation. That is, the com- 
puted procurement inventory level for an item does not provide sufficient carcasses to 
allow repairs at the computed repair inventory level. To solve this problem, Naval 
Supply System Command (NAVSUP) has made some changes. Under the model inte- 


gration, there is only one RISK formula: 
RISK = IC,D[ IC4,D + ¿FE 


where: 
C, = (BIDXC) * (1 - Xo 


Also, rather than using a procurement leadtime or a depot level re- 
pair turnaround time, it 1$ uses an average acquisition time as the horizon for computing 


the safety level. The average acquisition time (L,) is defined by: 
L, = (1 — BID)L + (BID)T 


It 1s clear that L, is the weighted average of the procurement leadtime 
and the repair turnaround time because D-B represents the quantity to be procured and 
B represents the quantity from the regeneration. 

‘6, Model Evaluation. For the UICP model, three major problems exist. 
The first problem 1s that LICP models is a continuous review inventory system. The 
Korean Air Force repairable management system applies a periodic review. The second 
problem is that LICP assumes that no interaction exists among failures of items. 
Finally, the optimal inventory level for UICP model is determined by minimizing annual 
variable cost (order costs plus handling costs plus shortage costs). However, the Korean 
Air Force does not have the same objective. 

e. Availability Centered Inventory Model 

(1, Background. ACIM was developed by CACI Inc. under the spon- 
sorship of the Ship Support Improvement Project, Naval Sea System Command, 
PMS-306, and aproved in March 1981 by the Chief of Naval Operations for use іп de- 
termining consumer level stockage quantities for selected equipment [Ref. 20: pp. 1-2]. 


However, this model was initially used as a part of a larger model (Logistics Support 
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Economic Evaluation Model) which provided necessary inputs and enabled comparisons 
with other Navy stockage policies. After this, successive refinements for the simplifi- 
cation of the solution procedure and associated computer programs have been made for 
several vears. 

This model was originally designed to calculate inventory levels for 
all items in the parts breakdown of an equipment and at all stockage facilities in a 
multi-echelon support system. Thus, ACIM is capable of computing levels for all ships, 
intermediate maintenance activities, and depots that use or support the equipment. 
However, this model is mainly used in the provisioning process to compute shipboard 
allowances to achieve specified weapon system readiness levels which have not been 
achieved with the standard protection level models. 

(2) Availablity Measure. ACIM recognizes that the purpose of a supply 
system is to provide sufficient support so that a weapon system is operational when it 
is needed. The terminology used to described this goal is operational availability (4,). 
ACIM defines A, by the following formula [Ref. 21: p. 5]: 


A, = up time / (up time + down time) 


= MTBF / (MTBF + MTTR + MSRT) 


A, may also be interpreted as the probability that the equipment 1s 1n an operable con- 
dition at a random point in time. Among the three factors of A,, the MTBF and MTTR 
are svstem parameters outside the control of the supplv svstem and are viewed as con- 
straints. MSRT 1s the only term which depends on stockage postures, it is therefore the 
one that the ACIM model focuses on to achieve a given value Ѓог А,. 


(3) Mathematical Assumptions 


1. All parts are organized in terms of equipment with a top-down breakdown that can 
be represented as an arborescent network. Any part may be totally consumable, 
totally repairable, or any mix thereof. [Ref. 21: pp. 10-11] 


2. Stockage and maintenance facilities are organized in a hierarchical structure ac- 
cording to supply/ maintenance flows which can be represented as an arborescent 
network as illustrated below. Each facilitv has a colocated maintenance and supply 
capability. The facility at the top of the structure is assumed to have an infinite 
supply of all items. 


3. External demands upon supply are stationary and compound-Poisson distributed. 


4. All stockage locations use a continuous review, (S-1,S) ordering policy. 


50 


Cn 


Mean Time to Repair (MTTR) is defined to include all equipment downtimes that 
are not supplv related. 


6. There is no lateral resupply (transshipment) among bases. 


Items repaired at anv location are assumed to be returned to collocated stocks for 
issue. 


—1 
° 


$. If. for a given facility in the network, the stocks are physically distributed in several 
places, it is assumed that the resupply time for direct customers (next lower ech- 
elon) is independent of the location. 


9. The order policy assumption precludes consideration of economies of scale for re- 
supply. That is, all ordering is on a one for one basis. 

(4) Model formulation. The goal of ACIM is to maximize the opera- 
uonal availability (4,) of a Weapon system subject to a given inventory budget. With the 
definition given above of 4, and with the view that MSRT is the only term which is af- 
fected bv the stockage decision, the developers of ACIM argue that the allocation which 
maximizes 4, 1s equivalent to the allocation which minimizes MSRT. ACIM therefore 
actuallv attempts to minimize MSRT subject to given constraints. 

Let 1 be an arbitrary item in equipment e (which may be e itself). Let 
ч = O represent an arbitrary facility in the support system and u= 1,2,3,...., represent fa- 
cilities at the next lower level, i.e, those facilities that submit items for repair directly to 


or obtain resupply from facility o. Then, the objective can be explicitly stated as follows: 


Find values for S,, which minimize D, for all user locations u: 


subject to pros = B 
kv 


a ere: 
C, = unit cost of item k 
B = given budget for spares procurement 
D 


S,, = the level of inventory for item k at location v 


= expected delay per demand upon inventory for for item i at location u 


іш 


D., is one of the components of M, which is defined as the mean time to return a failed 


unit of item 1 at location u to a serviceable condition. M, is given by: 


Ma c Dy c a 
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where 
D,, = expected delay per demand upon inventory for for item 1 at location u 


tu 


I, — mean time to repair item 1 at user location u 


іи 


D, in turn 1s given by: 


D, (LIE s pn E 


Хэ 237 
"Where: 
S, = stock level of item 1 at location u 
À, ^ expected number of demands upon inventory for item 1 at location u 


P(X:A,T,) = probability of X units of stock reduction for item 1 at location u 


T, = mean resupply time for item i at location u 
I, is given by: 


T, 9 (Ly * Ly) *(1 — Вы + Вы) 


m 
| 


, = average resupply leadtime assuming stock 1s 1s available 


i. = additional resupply leadtime due to expected shortage at the resupply source 


£ 
R,, = average shop repair cycle assuming availability of spares for items within 1 
R 


4, = additional shop repair cycle due to expected shortage of spares items within i 


Since ACIM assumes that MTBF and MTTR are independent of the 
stockage policy, minimizing D,„ is assumed to maximize 4,, (operational availability of 
a equipment at location u). 

(5) Solution Technique. Assume that initial values for S, are given for 
all items and locations. These may all be zero or some minimum value given bv policv 
or current assets. Then compute the MSRT's for all items. 

ACIM calculates what the new MSRT would be if one additional 
unit of stock were placed against that particular item. Subtracting the new MSRT from 
the old MSRT and then dividing by the unit cost of the item provides the model with a 
value which is multiplied by the item’s demand to determine a selection-rank. After the 


selection number is calculated for each item, one unit of stock is added to the candidate 


with the highest selection rank number. The operational availability (A4,) 18 calculated 
and if the target 4, has not been achieved, the model continues the same procedure. 

(6) Model Evaluation. As with the UICP model, three major violations 
of the basic assumptions of the ACIM model exist. The first is that ACIM assumes that 
vou have a continuous review system. The second is that ACIM assumes there is no 
lateral resupply. Finally, ACIM assumes that the stocks are physically distributed at 
several locations. Obviously, the Korean Air Force supply svstem does not require se- 
veral physical locations, even though the svstem objective is to maximize operational 
availability. 

D. SUMMARY 
Throughout Section B and C of this chapter, we have reviewed the various inven- 
tory models and theories for repairable items which are currently in use in the military 


services, or in the process of being evaluated. All of the assumptions and the math- 


ematical formulations of the models were stated and discussed. 
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IV. MODEL PROPOSAL FOR THE KOREAN AIR FORCE 
A. INTRODUCTION 


In Chapter II, we summarized the current Korean Air Force repairables manage- 
ment system and repair process. The general structure of repairables management in the 
Korean Air Force is similar to the structure assumed by several of the models described 
in Chapter III. 


In this Chapter, we describe a proposed inventory control process for the Korean 


Air Force. This process is based partly on queueing models, and partly the Wilson- 


Harris economuc order quantity model. 


B. THEORETICAL BASE FOR THE MODEL PROPOSAL 
1. Notation 


c = number of identical repair channels 
ц, = average repair rate per repair channel 


о, = average failure rate of an individual unit of a repairable item 


^; = random variable describing the steady 
state number of NRFI units of item 1 


P = steady state probability that there are n NRFI units 


2. The M/M/1/K/K Queueing System 

This model 15 а limited source queueing model in which there are only K cus- 
tomers. It 1s variously called the machine repair model, the machine interference model, 
or the cyclic queueing model. It is one of the most useful of all queueing theory models. 
One way to view this model is shown in Figure 6. [Ref. 22: pp. 186-190] 

Ihe population of potential customers for this queueing system consists of K 
identical devices, each of which has an operating time of O time units between break- 
downs, O having an exponential distribution with average value 1/a. Тһе опе 
repairman repairs the machines at an exponential rate with an average repair time of 
1 / time units. The operating machines are outside the queueing system (outlined by 
the dashed lines) and enter the system only when they break down. The queueing system 
always reaches a steady state because there can be no more than K customers in the 


system (one machine being repaired and K-1 waiting for repairs). 
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Figure 6. Оцнсисшгс Svstem wit; One Repairman (M/M/I/kK/[KX) 


When n of the machines are down (not operating), then K-n of them are oper- 
ating and the ume until the next machine breaks down is the minimum of K-n identical 
exponental distributions with parameter (K-n) o. 

3. Tne M/M/c/K/K/ Queueing System 
This qucucing system is similar to the machine repair model considered in the 


COME DOD 


section above except that we have c rather than one repairman, where c 


Мат” REI 22? pp Poorer) 


(C БЕЗЕГОРИАСА МӨРЕТ FORJEHE KOREAN AIR FORCE 
|l. Application Conditions 


ОИ 1С Пс ОО өтемесе available at all times. fn this case, 
ее с ег Olana luch the Korean ^ir Force wants to have 
айдаш Ойым ПЛ (рар; 


tJ 


P(Successful Repair) is close to 1. When P(successful Repair)= 1, we have a pure 
machine repair queueing system. When P(Successful Repair) 2 0. we have a con- 
sumable item Whose management is best accomplished using standard methodol- 
ogies [rom inventory theory. The sample data collected from the Korean Ai Force 
for this study indicated that P(Successful Repair) 2 0.98 or higher for all repairable 
items examined. l'or this reason, the proposed model i5 based partlv on the nxi- 
chine repair qucucing model. The data base used for the study was obtained from 
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Figure 7. Queueing System with e Repairman (M/M]Jce/K/K) 


the KAFLC during the 1984-1988. Vabic 3 shows the mean probability. of sie 
cessful repair based on the sainple data.6 


TiC — ICCESSFULR = COMPUTED FROM SAMPLE DATA 


Average 


1.000 1.000 1000 


0.980 0.987 0,954 
0.986 0.976 0.987 ШЕ 





бэ 


Conditions close to the SUNT RIS gueueing sesicn NOL но її 
number of repairmen in repair depot, 1t might be reasonable at first to assume that 
the MM cC KAR model applies. | llossever,opagstureparrsesequire «рес исо х 
held bv few people for a given repairable part. Further, the repair depot is stalled 
and scheduled so that there 15 litthe overall excess capacity. лба cs Tr ye n 
little change in the repair process output rate when the second, third, fourth, etc. 
unit backs up in the repair queue. 


dj. Steady state has been achieved. 


6 The sample data list appears in Appendix A and B 


сл 
о. 


2. Mlodel Formulation 

This proposed model consists of 3 steps. The first step computes a minimum 
population size A, for the ith type of repairable item. Population refers in this case to 
the installed material plus the stored RFI material plus the NRFI materia! awaiting re- 
pair. In other words, all material, both good and broken. 

Once the nunimum population size K, is chosen, an attrition buy safety level 
5, 1$ chosen. This safety level 1s used to insure that the actual population size usually 
stavs above A, during the period when the system is awaiting the arrival of an attrition 
buy. 

The last step is to compute the size of an attrition buy itself. This Q; is com- 
puted using the Wilson-Harris economic order quantitv model. 

In step 1, described on the next page. the МАМУ КК queueing system will 
be used to obtain probability expressions of the form PANT > K, — D). These expressions 


are obtained as follows: 


K, 
р> к-р) = 2, P(N-m 
nz K; x D, +1 
where: 
N, 7 Discrete random variable representing the number 


of units of item i which are broken 


K, = Decision variable representing planned minimum 
number of units of the item to be owned, including those 
installed in the repairable equipment 


D = Number of units of item 1 which the Korean Air 
Force wants to have available at all times 


The right hand side of the last expression is equal to: 


K, — D, 


] — 2 P(N, = n) 


n=0 


using the М/М/,К/К queue, Allen [Ref. 22: p. 192] gives the following formula for 
P(N =n): 
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where: 


К, 
Кү ши 
k=0 
In step 2, described in*the following section, tbe expression we ES IS alo 


used. It 1s obtained using the Poisson distribution: 


К “(0,7 
(Dail ye Pa 


P(X,= K) = a 


where: 


X, = discrete random variable representing attrition 
demand for item 1 during a procurement leadtime 


D,, = forecasted attrition demand for item 1 
L, = forecasted procurement leadtime for item i 
We will want S, such that P(.Y, » S) is less than a risk value, b,, which is pre-specified 


bv the user, where Р(Х > S), is given by: 


ОО 


PX»-s)- D P(X; = k) 


= (Dit) ee ae 


k! 
kz S; +1 
2 k —(D,L) 
21-7 Gute eh 
ОИ k! 
k=0 


a. Step I 
There are several different objective functions which might be used to help 
choose the size of the population fot the ith item. The following four subsections de- 
scribe four different math programs for choosing K,. 
(1) Alternative I 


Choose K, by the following: 
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min K 
subject to P(N; > (K; — D;)) < В; 


where: 


m 
II 


‚ = Number of units of item 1 which the Korean Air 
Force wants to have available at all times 


b, — Parameter chosen by the decision maker 


This first alternative formulation contains no explicit budget constraint. However, 
mininuzing K, for all items is equivalent to minimizing the cost of that portion of the 
population which 1s owned to meet the availability goals. Unfortunately, this formu- 
lation doesn't identify solutions which exceed the available budget. The next math 
programming formulation accomplishes this. 


(2, Alternative II 


nm 
min УХ, 
i=1 


subject to ANS Raz DISD cim ит 


m m 
and У ВК, + > айтах(0, К, — К.) < В 
i=] i=] 
where: 
h, = annual holding cost per unit held 
B = total budget for one year 
m = number of repairable items being considered 
K,, = previous planned minimum number of units owned 


a, = unit cost of item 1 


/3, Alternative III 


Maximize customer service subject to investment constraint: 


2 


m 
min Z = 2 P(N,> (K — D) 
i=] 


m 


m 
subject to У Ак, zi У а|тах(, K; — Кл) < В 


i=1 LJ 


(4) Alternative IV 


Maximize customer service subject to investment constraint (another form): 


min 4 = max {P(.\,;>(K; — D,))} 


n nm 


subject to » AK;  ? almax(0, K,— K,)] « В 


i=] i=] 


The choice between these alternate formulations depends on the way the Korean Air 
Force feels most comfortable with resolving the conflict between budget limitations and 
customer service. For the numerical computation of this study, we chose alternative I. 


This choice balanced the need for realism against the time constraints imposed on this 


research 
b. Step 2 
Attrition buy reorder point, R, = K, + S, where S, is obtained from the 
following: 
min 59, 
subject to P(X; > S;) < ba; 
where: 
X, = discrete random variable representing attrition 


demand for item 1 during a procurement leadtime 


b, — parameter chosen by the decision maker 
for item 1 


Using the reorder point R,, place an attrition buy when the total number of unattrited 


units (working and failed, installed and uninstalled) reaches or drops below А. 
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c. Step 3 


Use the Wilson-Harris EOQ model to choose procurement quantity: 


© NV IG 


A, = administration and ordering cost for item i 


where: 


D., = attrition demand rate for item i 
I, = annual holding cost rate for item i 


C. = unit cost for item i 


I 


So the initial buy quantity formula is R, + Q = K + 5, + О. 

This proposed repairable item inventorv control process works as follows. 
To start the system, purchase K, + S, + Qi units. As installed items fails, repair them 
as soon as possible (i.e. do not use batch repair). When attritions in the repair process 


reduce the actual population size to K; + S,, place an attrition buy. 


D. SUMMARY 

In this chapter, we developed a new inventory model based on the M, M/'I/K/K 
queueing theory and the current repairables management system of the Korean Air 
Force. In the next chapter, we will compare the current inventory model with the pro- 


posed model, in terms of the numerical results of each, using computer computations. 
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V. COMPARING THE COMPUTATION RESULTS OF THE MODEL 


A. INTRODUCTION 

In this chapter, we will discuss a specific result for the two models, one 15 a existing 
model of the Korean Air Force, the other is a proposed model. The results show only 
summary of the computer computation and a comparison of the results. The whole 
computation results are shown in Appendix C. 

The purpose of this chapter 1s to provide a numerical analvsis between the proposed 
model and the current Korean Air Force model. The data used in the computer com- 
putation consists of a selection of 8 line items form the data on 50 line items which was 
obtained from the Korean Air Force. Due to the limited time available for this research, 
only Š items were selected for use in the comparison of the models. A selection of high 
demand and low demand items were chosen. The complete data set is described in Sec- 


tion C of this chapter. 


B. MEASURE OF EFFECTIVENESS (MOE) 

The following MOEs are all reasonable methodologies for evaluating repairable item 
inventory models for the Korean Air Force. However, we will concentrate on one MOE 
for model comparison such as operational availability. 

l. Fill Rate 
It is computed by taking the total number of units demanded at a base over a 
fixed period time and dividing that number into the total number of units issued at the 
time they were are demanded. Thus, fill rate is the percentage of demands that are filled 
immediately. 

2. Operational Rate 

Operational rate is the probability that, at any given point in time, there will be 
no stockouts from base supply (backorders). Operational rate is computed by counting 
up the length of time (in days) during a year that no backorders existed and dividing this 
number by 365. This gives us the percentage of time during the vear that no backorders 
were in existence. 

Operational rate has an advantage over both the fill rate and mean number of 
backorders in that it may be directly related to the supply system's effect on operations. 
However, it has a disadvantage over both the fill rate and mean backorders in that it has 


a rather bothersome all-or-nothing character. 
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3. Mean Supply Response Time (MSRT) 

МКТ 15 the mean time it takes for the supply svstem to respond to the demand 
for a replacement part or component. It is obtained by taking each stockout that is es- 
tablished during a period of time, observing how many days it takes to satisfv the 
backorders, adding up all these numbers and dividing the sum bv the total number of 
demands during that period. It 1s the average time weighted unit short (TWUS). 

MSRT 1s important by itself as an indicator of success of the supply system in 
meeting response time goals. As a measure, it considers both the likelihood of satisfying 
demands from stock on hand and the length of the delay in satisfving demands when the 
svstem runs out of stock. 

4. The Average Number of NORS Aircraft 

NORS stands for “not operationally ready because of supply”. The average 
number of NORS aircraft is computed by counting for each aircraft the number of 
NORS davs during the course of a vear, adding these numbers up for each aircraft, and 
finally dividing by 365. Considering that the purpose of a spare parts supply svstem for 
the Korean Air Force 1s to maintain the operational readiness of aircraft, the average 
number of NORS aircraft would certainly seem to be a reasonable measure of effective- 
ness. A stockage model attempting to optimize with respect to a NORS measure would 
require more restrictive assumptions than those optimizing with respect to fill rate. the 
average backorders, MSRT. or the operational rate. Above all, its lack of “separability” 
IS à cause of major mathematical problems. 

Something to note is that the operational sector of the KAF has consistently 
suggested that the Korean Air Force supply system should be managed in operational 
terms, not bv supply terms. That is, their imminent question has been "How many 
spares are needed to provide an operational readiness of X percent of the aircraft?". 

5. Constraints 

The budget available for the investment in spares 1s either a constraint in all of 
the models or the objective is to achieve specified performance at the minimum budget. 
The essence of any inventory control problem is the trade-off between cost and system 
performance. As mentioned in the beginning of Chapter II, the budget constraint is the 


major resource constraint for the Korean Air Force. 
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C. COMPARING THE MODEL RESULTS 
1. Data Characteristic 

The data used for the model comparison are 8 line items of the 50 line items in 
the sample data. The whole sample data is shown in Appendices A and B. This data 
consists of: a) quarterly demand data for 19 quarters, and b) quarterly attrition data for 
19 quarters for each item. The demand data was analyzed using Minitab. The mean, 
standard deviation, variance and squared coefficient of variation were calculated for each 
item. Since data concerning the failure process for each installed item wasn't available, 
we estimated the mean item failure rate, «, by computing mean demand over 5 vears, 
divided by number of aircraft. We estimated the annual item repair rate, џи, by dividing 
actual working hours per vear by the average repair time in hours for the item. The 
current average repair time for each of the 8 items in the sample was obtained from the 
Korean Air Force [Ref. 23]. These values are shown in the table below. For the pur- 
poses of this comparison we assumed that the random variables associated with the 


failure process and the repair process had exponential distributions. 


Table 4. MEAN REPAIR TIME: 
This table shows the mean 
active repair times for the 8 
sample items 


Mean Repair 

Time (hours) 
|_| tna эт 
Гог 
TS TG 


22 | a 
51 аа 2s а 


и 





The actual data for an item consists of 19 quarterly demand observations during 
1984-1988. Items were selected from three different kinds of aircraft, the F-4D/E, 
F-5A/B and F-16A/B. Appendix B shows the Minitab results for all 50 items in the 
sample data set. 
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The data used in this research also included quarterlv attrition quantities for 19 quarters. 
However. for the model comparison, we selected 8 items.? The items which were chosen 
had a varietv of annual demands. 

2. Comparing the Results 

Ihe complete computation results for both models are shown in Appendix C. 
A summary of these results is shown on the following pages. For the Korean Air Force 
model, the target operational availability is 90 percent. In the Korean Air Force, this 
is a fixed level of availability. 

During the computational steps, the actual operational availability allows an 
increase from 90 percent to 93 percent, however the proposed model didn’t meet the 
budget constraint for a 94 percent operational availability. Figure 8 shows the difference 
of Kk values: a comparison of the current K values of the Korean Air Force with min K 
from the proposed model. Figure 9 shows the difference in safety levels for the proposed 
model and the Korean Air Force model. In the comparison, the attrition protection 
levels of 0.98 1s fixed. Appendix C shows the results based on several attrition protection 
levels between 0.90 and 0.98. SLQ is a safety level for the Korean Air Force model, min 
S is the attrition safety level for the proposed model. Figure 10 shows the difference in 
procurement quantities. For the Korean Air Force model, RO is a procurement level, 
optimal Q 1s for the proposed model. In the actual difference for the two models, the 
proposed model has slightly lower procurement quantities than the Korean Air Force 
model. 

3. Description of the Computer Program 

For the Korean Air Force repairable item inventory model, we first computed 
a demand forecast using a trend value to guide our choice of alpha value to be used in 
an exponential smoothing forecast for the next quarter's demand. When the trend value 
fell between 0.9 and 1.1, we chose an alpha value of 0.2, otherwise we used an alpha of 
anu ren tijewemand forecast Was Used to calculate the DDR, DRP,RC@, OSTQ, SLQ 
and RO. The computational procedure for this current Korean Air Force model 1s 
briefly described in the following paragraph. The complete computational results for the 


Korean Air Force model are shown in Figures 20 to 27, and the whole mathematical 


7 ]tem numbers for the 8 items chosen from the 50 items are 10, 16, 20, 27, 30, 33, 37 and 
46. 


8 Attrition protection level is the probability that attrition will not reduce the population of 
item 1 below XK, when an attntion reorder is placed when the population drops to the K, + S, level. 
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Difference of Optimal Procurement Quantity 
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process for the Korean Air Force repairable item inventory model is described in Chap- 
ег 

For the DDR, we computed the total demand for each quarter, divided by 365 
ux ine DRR equals the RTS divided by the sum of RTS, NRTS and COND. Ihe 
RCQ was computed by multiplying DDR times DRP times RCT. This RCT refers to 
the repair cycle time allowed for depot level maintenance. For these computations, we 
used 120 davs for RCT because this is the maximum allowed by the Korean Air Force. 
EOSTO is computed by multiplving DDR times NDRP times OST. For OST, we 
used 365 days which is the maximum allowed by the Korean Air Force. The SEQ is 
computed bv finding the square root of 3 times RCQ plus OSTQ. The RO 1s computed 
bv adding SLQ, RCQ and OSTQ. 

For the model comparison, we used two system parameters (SLQ and RO) and 
the item's current K value. Using the historical demand data obtain from the Korean 
Air Force, and their current repairable item inventory as described above, the requisition 
objective (RO) and safetv level quantity (SLQ) were computed for each calendar quarter 
in 1988 (note that the equivalent historical information was not available from the 
Korean Air Force). The current K value in Figure S was obtained by taking the sum 
of RFI, NRFI and the currently installed population, as of November 19SS (previous 
historical data was unavailable). The SLQ for model comparison in Figure 9 is the av- 
erage safety level in 1955. Also, the RO in Figure 10 is the average RO level in 1955. 

For the proposed model, we first need to choose 6,, for the first constraint set 
for each item. For b one minus the current target operational availability in the 
Korean Air Force was used. To obtain the constraining total budget value B, we com- 
puted the sum of the current K,, value for each item times its unit cost plus the holding 
cost for each item, over all 8 items. The rho value (+) for an item was derived from 
the average failure rate (a) of that repairable item, divided by its average repair гате (и). 

For the model comparison, a value for the budget constraint had to be obtained. 
We estimated a value for this variable by summing the following products over all items: 
K,, times unit cost plus annual holding cost for item 1 (reminder: K,, is the number of 
units of item 1 currently owned bv the Korean Air Force). 

In step l, the min К, values are computed based on two constraints. First, we 
calculate the probability that .V, is greater than K — D, for an item. The starting value 
used for K, is D, (the size of the currently installed population). For the min K, in the 
computer program. we first compute P, of the M/M/1/K/K queueing svstem. And then, 


ше compute P(IW— n) using Pa In the P(N, = n) computation, equals- K, —2D;. If*the 
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constraint on P(N, > K, — D) isnt satisfied by choosing K = РО шеп h re s ip 
1. These steps are continued. until P(N. > K, — D) is less than 0.10. At this time, we 
took the current n value, and then added D, to get the min A, value. After obtaining the 
min K, values from the first constraint set for each of the 8 items, the budget constraint 
is checked. If these K, values don't satisfy the budget constraint, we go back to the first 
constraint set and increase the b, values. In this case, a method for choosing an in- 
creased set of &,, 's must be chosen. Two wavs exist to do this. The first method 1s to 
increase the 6,, values simultaneously for all values. The other wav of solving this 
problem is to apply different operational availabilities for each item. In other words, the 
operational availability chosen would depend on the item characteristics. The first 
method described above (simultaneous increase of every 5,) was used in this work. If 
the min K, satisfies the budget constraint, this 1s min. АК, used the model comparison, in 
Figure $. As it turns out, the first set of nun K, values met the budget constraint. So, 
the min K, values from the first constraint set gave a feasible solution in the budget 
constraint. Actually, the budget constraint allowed up to 93 percent operational avail- 
ability for all $ items. 

For step 2, we use the Poisson probability distribution to get minimum attrition 
safety levels. Suppose that risk 6,,= 0.1. This would mean that the attrition protection 
level is 0.90. Therefore, 1 — 6,, equals the protection level,” For the min Sy cuirst cal 
culate probability of X,=0, and next calculate probability X,=0 and so on. We sum 
these individual probabilities up until the cumulative value meets or exceeds the pro- 
tection level (for the model comparison, we use 0.98). This is the min S, value for the 
model comparison in Figure 9. 

For step 3, the Wilson-Harris EOQ formula is used. For this formula, a per 
order ordering cost of 10 percent of the unit cost was used. This is the value which the 
Korean Air Force currently believes to be in effect [Ref. 1: p. 65]. Similarly, a holding 
cost rate of 10 percent per year of unit cost was used. The resulting optimal procure- 
ment quantities for the model comparison are shown Figure 10. The complete compu- 
tational results for the proposed model are shown in Figures 28 to 35. 

4. Summary of the Results 

The proposed model allows for a 93 percent operational availability within the 
budget constraint. This 1s an improvement of 3 percent in operational availability with 
the same budget that was used with the current Korean Air Force model to obtain 90 


percent operational availability. 
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D. SUMMARY 


In this chapter. we compared two models based on computer computations using a 
sample data set. The proposed model has shown an improvement over the current 
model. However, it doesn't mean the proposed model is better for the Korean Air Force. 
A tremendous effort would be required to actually implement this model in the Korean 
Air Force. Furthermore, additional research should be done to explore the performance 
of the proposed model. 
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VI. SUMMARY AND CONCLUSION 


For this research, there were three major objectives. The first was to review the lit- 
erature to understand the current Korean Air Force model and other models, such as 
METRIC, Mod-METRIC, for controlling the stockage decisions in multi-item, multi- 
echelon inventory systems involving repairable items. The second objective was to de- 
velop a new model for the Korean Air Force inventory management system. The last 
objective was to evaluate this new model by comparing the current Korean Air Force 
model with this proposed model. 

For the first objective, we studied the METRIC family model, ACIM model and 
L.S Navy model. We then reviewed their basic assumptions and mathematical ap- 
proach. For the second objective, we reviewed the Korean Air Force inventory man- 
agement system concept and conditions, and developed a proposed model. In order to 
accomplish the last objective, sample data from the KAF was used in comparing the 
Korean Air Force model and proposed model. 

Analysis of the sample data revealed that the proposed model vielded some im- 
provement in operational availability and safety level with the same budget. As men- 
tioned in Chapter V, the proposed model requires further research. Currently the 
research is hard to do because of time constraints and data available from the Korean 
Air Force. So, additional research should be done to eliminate these shortcomings in the 
proposed model. Also, for the implementation, we need more data points to obtain 


better distributions for the demand process and repair time process. 
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АРА. 5АМ Р ЕРАТАЫ&Т. 1 


A. SAMPLE DATA DEMAND DISTRIBUTION LIST 


Table 5. DEMAND DISTRIBUTION LIST 
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Table 6. DEMAND DISTRIBUTION LIST 
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Table 7. DEMAND DISTRIBUTION LIST 
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Table 8. DEMAND DISTRIBUTION LIST 
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Table 9. DEMAND DISTRIBUTION LIST 
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APPENDIN B. SAMPLE DATA LIST 2 


A. SAMPLE DATA DISPOSAL RATE DISTRIBUTION LIST 


Table 10. DISPOSAL RATE DISTRIBUTION LIST 
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мај ој ој ој ој ој ој ој ој ој о 
оош рни 
|. 883! 0| 0 


P(successtu] 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 0.92 
repair) 
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lable 11. DISPOSAL RATE DISTRIBUTION LIST 


отк пј mc, npj| nj i| 16| пј вј 19| 20 
Шаа ом] ооо 2 
КЕ 90395901590 0| 0|.2| 0| 01:01 0 
мај ој ој ој ој ој ој ој ој ој о 
Ж?» 9.01 9.5910 50) 0.0.21) 0 
54) 0) 0) 0| 0) 0) 0) 0) 0) 0) о 
082) 0) 0) 0| 0) 0) 0) 0) 0) 0) 0 
| ај ој ој ој ој ој ој ој ој ој ој 
ШЕ e ој ој ој тј ој ој ој ој of o 
— кор 
| 862ј| ој ој ој ој ој ој ој ој ој о 
еј ој ој ој ој ој ој ој ој o 
SS веај ој ој ој ој ој ој ој ој ој о 
ры, 90.9 0 09 @| o| oj oj o 
82) 0) 0) 0) 0) 0) 0) 0) 0| ој о 
| 823] oJ o[ oj oj o[| ој ој ој ој о 
874) ој ој o[ 1| 0| Ol oj oj oj t 
5 ој ој ој | ој ој ој ој ој ој 
— Ја 
f É šJ] 0) 0) 0) 0) о 


P(successfu] 1.00 | 1.00 0.92 1.00 
repair) 
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Table 12. DISPOSAL RATE DISTRIBUTION LIST 


| QTR | 21)! 22| 23| 24| 25) 26^ 27| 28| 29 | 30 
| ај ој лк 20 a 
| мај ој ој ој ој 1| 0| 0 то а 
ваз ој ој ој ој ој ој ој ој ој o 
|. 84] o] o| 0| 0|] ol ua MEM 
— 81| o Too о 20| 11| о 
8521 0] 01 0) 3) ој ој | ој о 
| 8533] o| o! 0| 0| Т| 1 770 ae 
[8| o[ oJ oj oy 0) о| oj o^ oj oj 
8611 ој ој ој ај ој ој ој ој ој о 

5362| ој ој ој ој ој ој ој ој ој 0. 
| 8531 0 ој 0| o| 9 5 xol ME 
_ 864 0| O 0| 2|] 509 лаға аа 
|. 8r1| oj o| oj o| mu өт” ч аа 
| | 812| 2] 9| oj o| 1| ој ој јони 
|. 8723! 0| 0j oj o, oj ој ој ој тј ој 
(ај ој ој 1| 3| ој 4 000 а 
| 1| ој ој ој чј ој 1 2 EOS 
لا‎ р 
|. 883! 0| ој 


P(successfu] 0.94 | 1.00 | 0.96 | 0.97 | 0.97 | 0.97 0.96 
repair) 























80 


Iable 13. DISPOSAL RATE DISTRIBUTION LIST 


ҰН а ы "ыы >> 





кз е о ое 


P(successfu 0.99 1.00 0.96 | 0.94 
repair) 
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Table 14. DISPOSAL RATE DISTRIBUTION LIST 


мар ог оо в 06 № M меси 
Га о of of of of of 1] of 1] 0 
Гэр е of of oF of of oof 















сєрцэрцэр эрэ сэрж 
sss] ol o 


P(successfu] 1.00 | 0.97 | 1.00 | 0.97 | 1.00 0.96 | 0.92 | 0.98 
repair) 
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B. SANIPLE DATA ANALYSIS 


Table 15. DEMAND DISTRIBUTION ANALYSIS LIST 


E | їй эн спој 286 | 0.79 
L5 | 19 | 3er | m | 3:9 | e» 
i | os 
рев ар чај гр ш- 
e | ue | 2900 Í ов | 
TN NE 





















СЛ 
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Table 16. DEMAND DISTRIBUTION ANALYSIS LIST 


2 | чај 21939 | — 369 | в зе 
so п | — 329 | —179 зе | ве 
H3 [ | омог 095 | 026 оа. 
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МАИГЕ C. ВОМ ОЕК PROGRAMS AND RESULTS 


program ROKAF (input,output); 
(уг? 

tnb10 өсгау [1..19,1..0] oF croal; 
tab 20 arroy (1. ES. x сов); 
tab In array (1. .19,1..0] of геа1; 
vnr 

clik * booleon; 

1.J).k : interer; 

in_tab : tabil; 

DOR Кар : (1210; 

DRC tab : tnb10; 

ПСО кә : tab)JO: 

0570. (лы : 18010, 

SIQ_tnb : 500910, 

ПО. Қа о: 

АШ" бар : tnb70: 

TRUD tnb : tnb20; 

ECST. tnb : tnb270; 

EXP tab : tnblu; 

DIS_tnb : tabio; 

TEWUP_tob : tabu; 


PROCEDURE titlie(tk:integer); 


begin 
caso k or 
l : wuritoln 
: ucltotn 
нгііліп 
urltcin 
uritcln 
uríitein 
writedn 
; writeln 
end I S и ) 
and; 


яғ” TREND ANALYSIS жат», 
11% EXCONEMEIAL  SHOOTHING 
trí DDR 10047 
18313 Di $*t*ttt" 
КЕРЕ нса даа 
wert 0514 кте” 
аа: 51.0 vere 
эн Ro 132447 


. . . . . . . з 


ЕКОСЕРИКЕ rec (l,j, ki integer;vor TENP_tob :tab3O ); 
begin 
caso k of 
TEMP tab[k, 
: TEHP_tnb[k. | 
* TEUP. tab(k,i 
; TEHC, tab(k,i 
ГОЛЫГ E des 
: TENE Neh ae 
ТЕП” (ГТО, 
TEUP_tab(k, 
us (case k) 


It! tnb [i 
EXP 
ООД tab 


пса tab (i 
0574. (а0( 


ци ни ии 


TM tab ШІ. 


3 0) -)) t5 & cC) м 


e 
end; 





!ieure I1. Computer programs for the KAF Model 


(****t*. MALI CROGRAI ttrt) 
begin 


(#>2>> DATA LECUT ROUTIIIE €€€«x*) 


uriteln (' DEMAND DATA 1! ^); 
for 1 := 1 (06 13 05 
begin 


urileln (” FOR Уи) 

[or j о: Ше do 

rend (in_tab(i,j]): 
end; 


weritelu (° DISPOSAL RATE !!! ^5; 
Ѓос its Vita to do 
begin 
uritelu (' FOR . DN ,i); 
for И 0 «ale 
read (D1S_tab [ i,J))); 
end; 


(+›››› ТЫЮ MIALYSIS <<<<+) 


бөлін 
for i += I to do 
benin 
TRHD. tnb (i] :z (CCIH, Lub( 10,i]*1I ,tab[(19,i1])*2)/CIHTHI, tab( 16,i1)!1H tab[(17,i]* 
Ill, tab( 10,i] 
+ ТИ Lab(19,1312)5; 
ir CRUD tab [i] >= 0.9) nnd (TRUD tobtil < 


= 1.1) Uhen 
АБР 02572113) 221122 


else 
ACE ena s DE 
end; 
ond; 
(*›››› EXPOURUTI AL 51001111110) 424442 


for i Zf] toll) do 
EXP tab [1.;i] := 11. 
for { := 2 to 19 do 
begin 
Гог j -= I to UO do 
EXP. tab [i,j] := CALP..tab £5] + сар Тт + (( 1 - ALP tab {j)) 


ж ЕХР tnb (11-1,321): 
тай; 


for İi := 1 Lo б) 45 


FCST_tab (iJ := (ALP. tnb (iJ * Ill; tob (19,i]) + ((1 - ALP_tab [i]) 
t EXP tnb (19,i])5 


(*›››› DAILY DENAND RATE €t) 
lor i :2 1 to 19 dp 
begin 


гос ә S= jJ. (e В Ge 


DDR.tab [ i1,J) :s: ІП tabir 072 91. 
end; 


Figure 12. Computer programs for the KAF Model 
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(*»»»» DEPROT REDCAIT PERCEI <<<<%) 


Ner 1 :- 1 to 19 до 


begin 
for g := [to COS 
if "BIS ESOSI. J) с UCU then DNE ғар (1,3) :: 1.0 
elr^ 
begin 
Г И Басири, 590 - 00 Mien DRIP tsb (1,5) :- 0.0 
els^ 
UM SBA) | +: 01520255 (1,5) / 19 286 11.31: 
enc; 
iC (DiS tab £1.53) = 1) №1 (10 Lab (1.3) = 1) ühen DORE tab (i,.j) 
end: 
(%>>>> UCA A OSTO SLU КО <<<«?%) 
forme i :- 1 to 19 945 
ber in 
Гог 1:: | oO ёо 
berin 
ИСО кар 11.2) :: Ш9 8 ре [i,3) х. ОКЕ ба (1,ј]- 120; 
USTU tab [i.i] :z DUll tnb (3,4) “ ( 1 - UnP tab (i,3)) * 355: 
а Gab [1).)]] :- сігі ( 3% (са tob |1,3) * OSTQzteb [i.3))0); 
RO. tab [1.1] := GLO_tab(i.g) + СО сав(1,3) > 0570. ғар (1,3): 
end; 
"ridi; 
5 › PRINT ROUTINE <<<«%) 


ені := Бекі “ек 

for k := 1 to A Jo 
borin 

urileln; 

uriteln; 

vrikoln; 


tille (k): 
ыгіігіз (° ОТК 1 2 
479: 


CODE a ----- 


Lot ) :- P vow do 
begin 
uurite ( “D i): 2 
fer 1 121710 2 408 
begin 
rec Ci. E Eli tob); 
write (^ ",TENP-t8b[lk,1,3)): 
eur; ‘ 
writeln: 
end; 


Figure 13. Computer programs for the KAF Model 
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yl ТО Lhen 
беліп 
uriteln; 
writeln; 
«rite (' ALP’); 
[or j := i to 4 do 
write (^ '",ALP tob (3)); 
writeln; 
end; 
if k = 2 Lhen 
borin 
нгішсіп; 
writeln; 
weite (^ ЕСГ’): 
bor a sz 1 to 7020 
write (° ,ЕС9ГТ20201:217: 
uwrite)n; 
end; 


urilelu; 

ига 111: 

нгіббіп (° QTR 5 6 7 
Gr); | 

мгЇ1с10 (° --- ее 


иге DOCE 
[or 3 :: 5 іо 0 «сіп 
berlin 
rece (j. jk ПИЧИ Lab): 
НАС O ЕИ“: 


cud; 
ыгішсін: 
end: 


ШІ. П-ШІ then 
bertin 
uritoln; 
4rileln; 
ritel’ ALP’); 
[ne d :: b bo D da 
> 
нгі еіп; 
end; 
re ие арно 
berfin 
writeln; 
‘ writedu; 
write (’ ECT 
Гог 2 :2 5 
write(" “ 
writeln:; 
ent; 


D do 


J 
to 
,FCST tabt3l]):; 


end; 
end. 


Figure 1d. Computer programs for the KAF Model 


program test; 
type 
array а = апау |). 8] of real: 
һгглу К = лггау [1..0) o€ integer; 
аггау_х = array [1..9] оГ resol; 
var 
Пегас ба, 
ӘГГІК 2 аг ЕК; 
mr РС: АОИ: 
grec: array oh: 
arr_rho : nrrəy_n; 
ОСГ ЮО srray_a; 
arce xX arpay x: 
arr_Da : nrcray_n: 
nrr Q : nrray.n; 
arr. 5 : nrray_K; 
J total jJ Ne irem no, X.a Lemp : integer; 
lact,sum, lU, exp. cosl.,prot level : real; 
do aPaljn : boolcan; 


canis. 
Budget = 19945000: » 
Ic" 0.15 : 
ПА. 


(t»5»» SIEP I <<<<%) 


ГипеКјоп get FOUCRc:ànterter; rho: veal):repl; 
Var ` 


k : іпіспеғг; 


беліп 
Кер сс 1; 
“әп °= Í: 
Бәс k се 1) Lo Ke «lo 
bertin 
Гог 1:= Кс demiko Ke-kil do 
wet -= [act ЛЖ Лю: 
| ТИШЕП (Г К=з .k, fact: | ,Гас!):) 
ит: = о бит k lach; 


(СС сп = `,,cCum);:) 
ЕСІ = |І; 
ewi; : 


ТО = Абит: 
uriteln( PO = еши): 
спа; 


Figure 15. Computer programs for the proposed Model 
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proccslure initjynl_ncvuvr nrrc_n:array.n); 

begin 
SC ЕИ 
are_al2] 
nrr_n[ 3] 
nrr_n[4] : 
arrolo] 
nrc_n[|G] : 


11670.93r1.1; 
2579. ОЗЕ E 
S354.U0r1.1; 
1100. АЕ 
210 соғ, 1: 
7307.05.1; 


Hou Hu! gà Hog an и 


sre lel | 1HD60 721. T: 
пес лЦОЈ = 510132 7111: 
еп; 


proccdiro display minK(arc KA. : accay FK); 


begin 
urito]n(" min Kl s ",nrr.KT1]); 
urjleln(* тіп K2 = art KIZ 
vurilLnlıı( тїп ЁЗ = “sarr KOI: 
uritolu( тілім = “.arc_K[24]): 
uritelu(^ и Ко = ‘п ВЫ): 
uriteln(" min AG = ° arc K(G]): 
urilcln( min 57 = ',nrr_K[7]); 
uriteln(’ mir KO = '.nrr_RK(0)); 

erl; 


function їг ронег(13, : renljpower : inteser):real; 
ber in 


if pover = D then 
ulh power со 1 
nlse 


11. power го PM, + nth power(DL,pouer-1); 
ewl; 


unction fac(mim : іпісиег)гінігрог; 
ber in 
1 num z 0) Шел 
је s= Í 
ojze 
luc ге ат f! Fac(mim-]); 
end; 


[immctiei get, SCIL: cenli;prol Jevel:ceal): iuLerov; 
vnr 

U som : топ]; 

i : interer: 


Беріп 
Раш. КАЕ 
1 0: 
repent ‘ 
бил := Копия в римог(1,,121(еур(-115) )/Гас(1): 
IL E 
unl.il P. 5um 
met S := iel; 
end; 


Е 
s 
›= prot. level; 


Figure 16. Computer programs for the proposed Model 
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Disewe- 17. 


procedure (tet, DL(Cvor orr DL : аггау_а); 


begin 
пт bi i) = ОВ 
вее_01.(2] := 0.052; 
ocr DLL]. ОИ: 
pcr DL(4] := 0.304; 
асе. 01.15] := 0.230; 
вее_1)1,(6[ := 0.056; 
атг 117] := 0.3217: 
вге_ 06[0) := 92510; 

end; 


proceduro disploy.S(arr,S : nrcəy_K); 


беріп 
ucitel "оос oil): 
оге ,ave_S(2l); 
иго arr aloy; 
urite(” ‘,огг_6[4)]); 
Grite( ` ace 5[5]); 
шиг сс: ALE OID: 
игИиМ (.  ,лсс 5171): 
uriteln(” °,oarr_S(8}); 

end; 


nroccJure rct. injlCvar nrr Rc,arr D 


begin 
age kell] сс 90; 
nrt TEL = ШО. 
orc 0207-2115 
пег 5 с(411-:5 194: 
nec. hcl5] :s FU 
nec. Kc[6] := 199; 
ПІГ ІСІ?!) *= ІМ: 
nere КОГО сс 120: 
nrc pez 0M 
npe DEM ss Rh 
DEEP = № 
erc Di4] :z UP; 
arr D[5] 22001): 
"mec DIDI = М: 
ие. 1)[7] := 175: 
ег ПОЈ = ВВ; 


nrcc_rho[1] : 
nrr..cliol2] : 
ncc. rho(3] 
arc rlioof4] 
nrr, cho[5] 
acr. clio(6] : 
acr. chio[7] 
nrrc_rcho[”) ] 


0.00069, 
0.01010; 
0.00661; 
0.002371); 
000004: 
0.00652; 
0,005290; 
0.00071; 


Computer programms for the proposed Model 
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: агтду_К; vor orr_rho : 


агтау_л); 


Figure 15. 


ВР [1] = ШУШ 
агг (21 зе (0291. 
arrzx[o] 22102328 
arr. x(4] :- 0.93; 
arr xlo :s 02348 
nrr.x[(6G] :- 0:955 
өгг.х1(7) : 0.96: 
arr_x[0] := 0.97; 
nrr_x(9]J := 0.90; 

end; 

begin 


іһібілі а(пгг а): 
get init(nrr, Kc,nrr,. D,arr. rlio); 


rcpeat 

do »9nin :- false; 
iiem no :z 0; 

exp enst := U.i 


repeat 
ilem no := ikem по + J; 
СО := get PU(arr Kc(Ciblom no],arr rho( item noJ); 
loct =a: 
SIM = 
11:2220)| 
repeat 
Гог 1:= arr_Kefilem_uo] downto arr_Ke(lilem_noj-nti do 
[act := Cock 3: j4ner_rho[ item n5]; 
[ас з= [nctil ü: {necessary to comment) 
write(' i-o sm): 
sum:- sum * facl; 
uriteln(' sum - ',sum); 
[пс 221: 
n: n. 
mil. 3] €] - cum) «2 () - OA); 
uvileln min K = `,arr DI item uo ]tin-0); 
ore Ki iten nef := arr Df item uo] + a =|; 
exp Leost := exp_ensl + агг n[ ilem nol! (arr Df iLen i011); 
unLil item no = 8; 
if exp cosh <= Budget. Lhen 
display _minklarr_K) 
olno 
do nrain :z2 true; 
until not do nraiu; 


(t»»»» SIU TI <<<<%ғ) 


uriteln(' prot. 1evc] S1 52 53 51 55 56 57 587); 
нгібсіп(” --------- 2 ---------.............. D: 
for x :z 1 to 9 do 

Беріп 


Computer programs for the proposed Model 


Орч (0 arr_x[x]); 
ites_no := 0; 
got.DLCarr DL); 
repost 
item_no := ilem_no 4 1; 
nrr_S[item_no] := tok. SCare DLLitem no). arr x(x]); 


until item_no = 8; 
display.S(orr. 5); 
end; 

(™>>> — SIEP III <<<) 
ATE USGI := 0.2: 
пг [2] := 0:2: 
BEE UNI зе U.2: 
arr БМ | := 0.4; 
пег iba] = 0.2: 
re Valo: = 0.2; 
ас2200171::5:0:6: 
әгі. ОӨл|0) := 0.6; 


uritc)n; 
uriko)n; 


uritcin (“корзо ОСТА, PROCUREMENT QUAUITITY  «««ct^); 
Гог > 1201 іш Jo 
berin 
пгг 0 [x) :* zact ((2* (0.1 * (arr. 0[x]/1.1))* avr. Da 
urito (^ С: 
writeln, = '`,ncr_Q[x]): 
eid: 
urilein; 
нее НЕ; 
мг) 1.2] 1 (° К 5 Q TOTAL : 
uvilLcln ( --- --- --- 2 ----- 207 
for x := 1 Lo Ü do 
berin 
jC arr Alx] <= 2 Lhen Q Lemp 222; 
il nre Alx] c= 1 Июи q temp га; 
tola} :2 arr. Kix] + arr lx] + q leap; 
urite(’ E 
il x= 1 bhen wrilel’ °); 
vril.c( ` nrr.k|x]):; 
urjtc(' “улст ЭН 
write” a temp); 
1 x = | theu urile C ); 
«cile( ^ “korale : 
uriteln; 
enl; 
сік]. 


2 . ( 
Ele 19. 
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(x))/41C * (arr. alx)/1.12)2); 


$4191 
QTR 
D! 
p 2 
D3 
04 
WS 
00 
07 
pn 
109 
010 
U11 
Di2 
1J) 3 
914 
D) 5 
DIG 
D17 
DID 
1113 


< G O O — > 


C) n9 C 


озса са С», 


ALP 


атп 
01 
02 
pa 
U^ 
5 
DG 
1(7 
ро 
DS 
010 
011 
912 
113 
14 
(1. 
1106 
117 
рап 
119 


СОС САД но — N — ОФ 


NNC = == 


— — 


ALT 


Пе) 


Т 


ДЕКО ANALYSIS 


а ж» == ©» ж» ж» «а ер «ж» «е» шь s. 


. 00000000001; 00 
UVUVUNUNUDE нат) 
.Q0000000U00F ! QU 
ОООО ЕСО 
. OUOQU(OQ9( JOE O0 
. Qttguotrcioticg + OU 
. 0000000001100 
.0000U0000U0U1E (UU 


- 
— 


(U1J00100U0U8 F ОО 


ОНИН O0 
. UUDTUOCOOUUDIZ TOO 
. UUOUOQUDOUCGCE » 00 
Олег ОЗ) 
, UCOUCUOIDOGOCE TOO 
‚оогон o0 
. QOCQOttOOUUtE 00 
ООО ЕСО) 
.UU(0000O0000O0018U1 01) 
 OI(00U000000UE#+ 00 


.0000000000£E-01 


Со аа ан == — ағ ас «е Se ғ» ж» ж» ш» 


„ 0Оридцодококоо 
. DUECUOUDOTGUOO0T 1 00 
.UU0010U90UUUF 1 00 
. UO duioonco X 00 
. CHO THO OODLE (СО 
ЧОООН КОО 
ООО (t 
лаци СО 
“00100100 OG 


. UOOUt coo s on 
,UUO0QUUGOUOQE!I 090 
.UUungesceed E 00 
.UüUttagtajo oc од 
. Quay Unam Y 00 
ТИЛИННИН QI 
QOD CINE FO 
имин Ea 
‚сиң к) 
ОТОН ТТТ 


4.00009t0)0nmtz-Q01 


ссу — N COC = — 


N `N — C >) S O O tS 


.00000U0000F1 00 
.GUD0CCOOUCOOE rg 
.00000U0UO0UE ' OU 
.UUUOUUECUOUOUE C 00 
.0u 00000000 (UU 
.00000U0000X t UU 
.0000000000 00 
.ÜUOOUOGOHOOUE 1 UO 
. Q0UOQUUOUOE "до 


OUOUCN—O SO 


, UOC(COOCOGQOOE 1 OO 
. 01100009000 
.0OUItOU000001E 100 
. OtelggaooQoQggE gu 
.UOUtIU00000E 100 
.000000000QE100 
.QUUODOOOQUE 1 00 
. (0000QO0UUPDF ЈО 
. UD UDCUOOQUE 1 00 
.000000t(000H5 +00 


2.0000000000E-01 


.0000000U0U(QE 100 
.0000000000Е100 
ОООО 00 
ОООО ОЕ 00 
ШЕТТЕГІ ОГ 00 
ОНОН «uto 
¿f(U00J01C00)U010 1 GO 
, Ucucpooonemdpds 1 00 
цопадцоцос + OO 


2.0000000000E 100 


, UC OUGOaOO0DIGE 1 00 
, Unrated à un) 
рачуна до 
ПОРНО цо 
UOMO ey 
OUI ЕО 
OOOO 1 (10) 
«аладан 8 CU 
0 . UOoc0camcaodE 1 UO 


2.0000Unp009)]79-0 1 


t^ > a O — G O = =— 


C N — N N OO 


0 


энэ = — аһ чо аһ ар аһ со со کت‎ еә hm 


.0000000000Ек00 
.0000000000F100 
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п 2 2 en =5 $.23255059' IC 17-04 п г 9 sum = 7.61025390150:-01 
n = 1 nam z 7.E0025936/3110-04 n * 10 snm = 0.0191220000Е-01 
пос 5 sua z  1.05547423502T0.-02 n = 11 сим = 0.37236411455-01 
1227 eun: —1.410151867057E-023 п = 12 sum 2: 0.6715202772Е-01 
а :/ чал 20124652081 0725-07 и = 13 eun : 0.920100572385 -01 
n = f епа = 2.4n534510G4E-07 ° n = 14 rum = 2.1410221195£-n1 
поч rus z 3,2U020222^02E-02 min K = 102 
n = 10 eun : 4.0031725224E-03 РО = 2.0106501591Е-05 
n = 11 sun z 5.1600657356Uu6(.-023 а:0 sum = 2 0106501591£.-05 
n = 12 “ємл 206.4711724220Е-032 n*1 sum * 06.5310035056Г-05 
"noc 12 гии = 0.0562222694EF.-0U23 nz? sun z 1.1399011247F-04 
n = 14 rum 5 Ч.9034440:704Е-0) 12:23 sum 2 1.7740195344Е-01 
n 5 15 eum c 101.22425059210E-02 nos sum z 2.5902418201[.-01 
п 210 єнтө-:001.4959725729Е-02 n: 5 sum 5 3.661323 10558-0141 
n : I" 6 тав 201.91716911901Е-02 n z Ü sum = 5.02G0D0426357F-01 
n = 13 eun = 2.1950225%04F.-02 n = 7 sun 5 6.77123011921£-01 
n = 13 лип = 2.GJ72230275F-02 п: 0 sun z 0.9007002793Е-04 
n z 20 nun 5 3.1515065002FE-02 n = 9 cum = 1.17924G0530F-03 
п: @ run z 2.7465611G007£-02 n = 10 сне = 1.5310965239Е-03 
л 2222 sun = 4.4209507СҒ2;-02 n * tl sun = 1.9231217301F-03 
73 зим = 5.2135131336E-02 n = 12 sun = 2.5222597211Е-03 
n = 21 enum = 6.1032239396К-02 nose сал г 23.2020550001:-03 
п = 25 sun 2 7.10011727600Е-07 7. в = 14 сап 2 4.0391075037Г.-иИ3 
n= 2 sun = 0.230020 Иб5НЕ- 02 г 15 сит = 5.06122249 0 LE -02 
n г 2 cue 2. 9.5012050502E-02 n * 16 саг 0.31222501498-01 
n = 20 sun z 1.0903117057[.-1U1 n z 17 sum z 7.0250049147Е-113 
n s 29 sun: 1.2150556472£E-01 n= 10 sun =5 9.040274117Л5-03 
п + 30 sun z 1.4140199U35E-01 n z 19 sun = 1,1032101-191:-02 
а По 11 sum z 1.6000261764£-01 nz ela sun = 1.4134101921ЁЕ-02 
nz: 232 enn = 1.0007854160E-01 n s 21 sum = 1.751)1700018-02 
n = 33 sum 5 2.01707067096-01 n z 22 зэв 2 2.1151003099Г-02 
n 3 ТІ sum 2 2.2406U90133E-01 пап 5 22 sum = 2.$109G117410-02 
n 3 25 aun x 2.4151727072€E-01 n z 24 sum = 9.0970410559E-02 
n s 20 zum 2 2.755053J0157F-01 n= 25 sun = 2.6117525207Е-02 


Figure 28. Computational Results for the proposed Nlodel (OA — 0.90) 
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5.5221600918Е-01 
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7.7998160093E-01 
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.0721002019Е-02. 
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n =“ 7 sum = 4.045762209ПЕ-01 
n = 6 sum = 4.9610226544Е-01 
n = 9 sun = 5.003070$233Е-01 
n : 10 sun z $5.5190215420E£-01 
n = 11 sum = 8.0152270595Е-01 
n z 12 sun = 0.4572748720Е-01 
n = 13 sum = 0.07203023017E-01 
n z 14 sum z 7,2606192071Е-01 
п = 15 sun = 7. 0194097009Е-О1 
n = 10 eun = 7.9450250062F-0O1 
n а 17 sum = 0.2420200170Е-01 
n = 10 sum 2 0.507944902СЕ-О! 
n = 19 sun = 6.744¢41800233£-01 
n = 20 sum z 8.95201451861Е-01 
n = 21 зим = 9.1330300691£-01 
nin K г 109 
өшіп КІ > 43 
min K2 = 115 
main KJ = 234 
nin K4 = 100 
nin K5 * 102 
nin KO = 239 
nin K? = 200 
nin КО = 109 
prot. level S| S2 S3 54 55 56 57 50 
9.0000000000E-01 1 0 1 1 1 0 1 2 
9.1000000000Е-01 1 0 1 1 1 0 1 2 
8.2000000000Е-01: 1 0 1] 1 1 0 ! 2 
9.3000000000£-01 101 1 101 2 
8.4000000000Е-01 1 0 1 } 1 Q I 2 
.5030000000Ё-01 1 1 1 2 1 1 1 2 
9.0000000000Е-01 1 l 1 2 | | 2 2 
9.7000000000Е-01 1 1 1 2 1 1 2 2 
9.0000000000E-01 ] 1 1 2 2 1 2 ə 
*»»»»  OPTIHAU PROCUREMENT QUANTITY «cece 
Qi = 6.3245553202Е-01 
42 =: 565.3245553203Е-01 
83: 6.2245552203Е-01 
ОД г 8.94427198100Е-01 
45 : 0.3245553203Е-01 
06 = 65.32455523202Ғ-01 
07 = 1.0954151150Е.,00 
49 : 1.0954451150Е+00 
К 5 9 TOTAL. 
1 8j 1 1 n$ 
2 115 l 1 117 
3 234 1 1 220 
4 100 2 1 103 
5 102 2 1 105 
6 239 1 1 241 
7 200 2 2 204 
6 109 2 2 114 
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Computational Results for the proposed Model (OA = 0.90) 


ооо ER‏ ا ل . Y‏ و ید ن 


PO: 1.5048167631Е-01 nc: M sun z 23.0290254127F-U1 
" 2.0 сие 2 2J.50409n76231£-01 n = 30 sun * 12.215950:027Е-01 
п "ug = 5.027/304772E6E-U1 n = 39 sum = 2.G61223)5177;0E-O]1 
тои cum = 7.011)006%92232E-01 а: 40 зим = 2.9191)07 27172-01 
n “ 3 smn = 0.322220 ЛОЛЕ - 01 n = I cun = :1..271912167Е-01 
n z 4 чле 2. 8,0527030326E-01 n z 42 sum z 4.55213923033F-01 
n z 5 sun = 9.2520974П00Е-01 п : 43 sum z 4.87%0012602С-01 

min F. z: nj n = 44 зим 2 5.1995627701)Г.-01 

0 = 1.5011756G016F-02 и = 45 sum £ 5.52J12779191F-01 
n = U cum 2  1.5041/560466-02 "u = 46 спе: 5.0415402091Е-01 
n = 1 лот 2 012.3904721962Г-02 n 2 47 зим = 86.1000004421К-01 
n s: 2 бол = 29.1500022575[L-02 n : 40 зим = 6.470191600123!.-U1 
н z J mun = —7,.7504165220€-042 n z 49 nun z 0.77001760 12E-UI 
no ا‎ cum = 1.0220229521Е-01 n = 30 им = 7.000929^40:f.-01 
n = 5 suu 2 — 1.2156906521£-01 а г 51 зил = 7,3390074516Е-0! 
п т б sum = |1.б24017/2200Е-01 ма 52 sun 2 7,60260222С02Е-01 
а г 7? sun = 1.95711G6UUJ45F-01 п т 52 sum z 7.0530629015€E-01 
neon sun = 2.3173000431ЁЕ-01 n z 54 rum z 0.0010051539Е-01 
n* 9 sum = 2.069723129706Ff.-UI n z 55 sun * 6.207514222096Е-01 
n z 10 сия 2 2.09620200292Е-01 n z 56 sum = 0,5100754527Е-01 
nos 11 sum = 23.5112372404€-01 n г 57 sum т 6.696159U019E-01 
n s 12 зим = 2J3.9306566005E-01 п = 50 sun = D.86590069230E-01 
i-e 1j sum z 4.2741005002г-01 пе ОО sum = 9.0196130777€-01 
п = 14 sum = 4.01352370105E-01 n z 60 gum z 9.1577955017Е-01 
nass sun = %.252101)9502С-01 mln K z= 2235 
n = If num г Ó5.6E055U02556F-01 ГО = 92.2017781102Е-01 
n = 17 sum = 6.10920202327Г-01 n z f) sum = 2.3017791102Г.-01 
n : 10 zum 2 6.5190470:90Е-01 n = 1 sum = 5.5274127АС0Е-01 
п =: 19 sum = 0.9111293237Е-01 пл? sun 5 7.01201701030Е-0) 
п = 20 eoa = 7T7.42U02125tt0656-01 n = J sum = 0.05110762021-01 
n = ? sum = 7.0205141217C-01 n = 4 sun = 8.7225004921L-01 
п 22 "un 2 1.95099239906-01 n z 5 sum = 9.1672969001Е-01 
Па aun ® 0.2452007430С-01 nin K = 100 
п > 21 сим = 0O0.511415170^2E-n1 ГО = 1.0364105001Ё-01 
п = 25 sum = 0.149$110199Е-01 n = QO sum = 1,03041050911:-01 
а = 2 cum 2 9559900117 739F 11 n = 1 зим = 2.0214159“50(-Ц1 
n : 2? nun 5 03.14311007724Е-01 ns 2 бит: 2,91402443711:-01 

rin K г 11“ n z 2 sun е  J.0116503077£-01 

PO: 0.00101714059С-05 n z 1 аит 20 4.01407741008-01 
n = U зил = 8B.U)1047"596-0*5 nace зим = 5.31/750232n5E-01 
n: I sum 2 2.116909 SE a= & nun z 6.0126271199-01 
noy пим 2 2.4595: 1] 02110E-04 u = f num = 6.Г109%719ЛГГ-01 
nes 3 гөз : — 5.975800 41-04 ncn sum = 7,.1414124539715510] 
n * 4 сим z T.65n0gmor Wt -14 и = 9 rum : 7.01025391150Е-01 
n > nun 2 — L.0654/42325025-0U3 n * 10 sum z D.0191220006^5-01 
n = Û en 3 1.1709010 707-97 п = (І sun = 0.2721Г1411452-П1 
п 217 eum = [.90553361025-03 n s [2 тиа = D.67155u272/n6C-01 
nag гаш 2 RA. МЛН Е- 02 n z 13 сөт = 0.930410)90JJUUF -3)1 
п 21 ене о: 23.2102022225U08E-U: а: 14 sun = 9,14102Л41425-01 
n = 10 eum = 4.10:1725224E-U32. min RK = 102 
n 91 sua = 5.165067)5"05E-0:1 PU z “2.D1U6501571E-05 
n = 1?! sum * 0 .1711724228Е-92 п 220 свпе 2 2,010551)1541Г-017 
ІІ сөө 200.09529130114Е-032 "21! sum = 6.523100295156[;-0n5 
п = 11 “um 3 Я. 96241167042 -03 N sum = 1.139904117170-04 
п * 15 ntum = 1,271215"93101-1(2 и тд sun 2 1.,774019051411-01 
п 210 ‘зим z 1.49597257211[5 -102 n = 4 sum * 2.5902116201E-t11 
n = 17 sun = 1.01216G9011911F-02 а: 5 sua 2 132.0611234159Г-01 
ne: t eun * 2.19500255U041:-02 n 6 сит т 95.02004261571-01 
n =5 19 nun * 20172709007-2 n т 7 sun = 0.7713011:21Г-01 
п 5 29 um z 2.15160С651022-012 п. 0 san > 9.00075527921 -01 
а : 21 “ат = 3.745n11607E-U2 п зим = 1.1792440%30:-03 
в = 12 sum 2 4.,42090507604:-02 п = 10 gna = 1.5J1096523)7E-02 
п 22 sim =5 5.21251212265-02 п 21! sum = 1.9731217311ЁЕ-03 
n : 21 сим = 60.]U2372223!29 65-02 n * 12 а г 2.57299 P IDIE S09 
n: 25 sum = 7.1000/290008-02-., n s iJ oun = 2.2020550091Г.-1)3 
nz ?W син т 0.739412070308-12 n z 14 пал 2 4.0U2191U2Z59022Ff.- 03 
а г 71 zum  09.501215061121-12 n z 15 su z 5 0017J3J3)123111F-U3 
n z 20 rum z 1.0207117067Г.-01 n = 1G sum 2 0,3122250149Е-03 
n se un = 1]1.2150556472365-Ul n z 17 sum z 7.02566549117/C-02J 
n z 30 sum = 1,414041%4025Е-01 n = 10 коп 2 9.56402741191Г-02 
n : Jl зим = 1.GUUUIJ61761f-01 и = 19 зим = 1.1032161590Г-02 
а: 32 гии = 1.00072541805E-D01 n = 20 zun 5 J.41311015210-02 
0 = 33 enm =: 2.U1I70706709E-UI a:a sun 3. 1.7210170604€ -02 
n z 34 sum = 2,240609U123E-01 n : 22 на = 2.1151002339Е_02 
n z 35 sum = 2,4951727072Е-01 п: 23 suo = 2.5109611741Е-02 
n z JG tun : 2.755051U157U-01 n * 21 зим = 2.03704405595Г-02 


ricure 30. Computational Results for the proposed Model (OA = 0.91) 
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77292220072Е-01 
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.722595n n 01 
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" 2 5 sun z 2.09500021207-91 
n 6 rum z 2.22221671321-01 
n ? сит = 4.045/П220У0Г.-01 
n = Ü sum z 4.5010276544Е-91 
n 201 зов = $.0823470522323E-01 
n * 10 sun = 5,5195245426С-01 
n * 11 sun = 0.01523666U5C-Ul 
nose sum = 6.49/2746720Е-01 
n = 13 чна = D.0723n909D1/E-0U1 
n * 14 sun = 37.2C0619230/1E-Ul 
п г 15 sun = 7.61010070101-01 
n = 16 sun :. 7,945525000?2c-01 
пе 17 sum z 0.2420200170ЁЕ-01 
n * 10 sun 2 0.507941902ГС-171 
122219 тиа г 6.74410бб73ЛЕ-0 
n = 20 sum г 0.9520145101Е-01 
n = 2! sun =: 9.1220200694Е-01 
nin К = 109 
min КІ: 93 
nin К2 = 115 
min KJ = 235 
тіп К4 т 100 
nin KS :102 
nin K8 = 239 
min K7 = 20| 
nin KO = 109 
prot loval 51 52 SJ S4 55 S0 657 50 
9.0000000000E-C1 1 0 1 ! 1 0 | 2 
98.1000000000Е-01 1 0 1 ! 1 Ü 1 2 
8.2000000000Г-01 1 0 1 ! 1 0 ! 2 
8.301:31)Й/000001:-01 ! 0 1 1 | 0 1 2 
8.40U! mmi000£ -01 1 Ü 1 1 1 0 | 2 
g.5ygguüuUnDdungt-01 . 1]! 1 12552. ТВ | 2 
9g.6g000UQD0(QUE-01. | 1 1 2 1 1 2 2 
9.7000000000E-01 1 1 1 2 1 1 2 2 
9.0000000000Е-01 1 1 1 2 2 1 2 3 
*»»»» 0111НАН ГПОСУПЕНСИТ QUANTITY ccecce 
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Figure 32. Computational Results for the proposed Model (OA = 0.92) 
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